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Abstract 


A series of adaptive-wall ventilated wind tunnel experiments was carried out to 
demonstrate the feasibility of using the side-wall pressure distribution as the flow 
variable for the assessment of compatibility with free-air conditions. Iterative and 
one-step convergence methods were applied using the streamwise velocity compo- 
nent, the side- wall pressure distribution and the normal velocity component in order 
to investigate their relative merits. The advantage of using the side-wall pressure 
as the flow variable is to reduce the data-taking time which is one of the major con- 
tributors to the total testing time. In ventilated adaptive-wall wind tunnel testing, 
side-wall pressure measurements require simple instrumentation as opposed to the 
Laser Doppler Velocimetry used to measure the velocity components. 

In ventilated adaptive-wall tunnel testing, influence coefficients are required to 
determine the pressure corrections in the plenum compartments. An influence co- 
efficient is defined as the relation between the change in the flow variable at a point 
at the control level and the change in the pressure in an active plenum compart- 
ment. In the present study, experiments were carried out to evaluate the influence 
coefficients from side-wall pressure distributions, and from streamwise and normal 
velocity distributions at two control levels. Velocity measurements were made using 
a two-component Laser Doppler Velocimeter (LDV) system. 

Side- wall, top- and bottom-wall pressure distributions, velocity components and 
plenum pressures were measured systematically to evaluate the effects of Mach num- 
ber; the presence of a nonlifting model; location of the plenum compartments and 
the control points; and suction and blowing on the influence coefficients. Linear- 
ity and superposition of the influence coefficients was investigated. Free stream 
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Mach number characteristics of the test section were explored. Influence coefficient 
matrices were formed for both suction and blowing. 

A one-step convergence scheme was developed and applied theoretically and 
experimentally by using the influence coefficient matrices. The one-step scheme 
reduces the testing time by eliminating the intermediate experimental iterations. 
The success of the one-step scheme depends on the accuracy of the influence coeffi- 
cients. It was shown that the number of iterations could be reduced to a minimum 
with the application of the one-step scheme. Results showed that the same un- 
confined conditions obtained with the streamwise and normal velocity components, 
could be achieved when the resultant velocity distribution from side-wall pressure 
measurements was used as the flow variable. 

The present study was conducted in the NASA Ames 23 cm x 11 cm indraft 
wind tunnel with a ventilated adaptive- wall test section. Experiments were carried 
out at velocities from M — 0.5 to M = 0.75 and at angles of attack from a = 0.° 
to a = 4.°. 
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Chapter 1 


Introduction 


The testing of aerodynamic models in wind tunnels has been the primary means of 
obtaining information for the design of aircraft for almost a hundred years, begin- 
ning with the early attempts by the Wright Brothers in the 1890’s. Despite this long 
and largely successful history, it has always been recognized that wind tunnel data 
must be analyzed and corrected before it can be applied to flight conditions. These 
corrections must take into account the influence of the wind tunnel walls and the 
model support devices which are not present in flight, in addition to the changes in 
viscous effects that result from the change in scale from tunnel conditions to flight 
conditions (the so-called Reynolds number effects). 

In recent years, with the introduction of the fast computer, it is possible to deter- 
mine the flow around simple configurations from purely theoretical considerations, 
particularly when the influence of Reynolds number is small. However, in most 
practical applications, aerodynamic design continues to depend on both computer- 
derived and wind-tunnel-derived information. For this reason, it is essential that 
accurate means of correcting for wall and model support effects be incorporated 
into wind tunnel test procedures. This is particularly true in the transonic regime 
where wall interference becomes more significant and the influence of viscous ef- 
fects, through such phenomena as shock-boundary layer interaction, become more 
difficult to describe in a purely theoretical model. 
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Over the years, rational schemes have been developed and used to correct the 
wind tunnel data^ 1-7 ). However these schemes are generally approximations, based 
on potential flow methods, external flow field estimations, wind tunnel idealization, 
etc. Although they represent valid mathematical solutions, they may not describe 
the flow in the wind tunnel in all its complexity. Inaccuracies in wall interference 
prediction arise from the following facts: i) the governing equations are nonlinear 
in transonic flow, ii) wall boundary conditions in the ventilated wind tunnels are 
complex and difficult to predict, iii) it may not be possible to model the wind 
tunnel geometry properly, and iv) the flow may deviate from the two-dimensional 
flow because of the boundary layer on the tunnel side-walls. 

Since accurate corrections are almost impossible in the transonic regime, the 
wind tunnel flow should approach, as nearly as possible, the free flight conditions 
around the model. Conventional slotted and perforated tunnels are intended to 
achieve this goal, but unfortunately they are often inadequate, especially when 
models with high blockage ratios have to be tested. 


1.1 The Adaptive- Wall Concept 


The adaptive-wall wind tunnel concept is an alternative approach to eliminating the 
wall interference or to reducing it significantly by controlling the flow near the walls. 
The principle of the adaptive- wall technique can be explained easily by referring to 
the flow field around a model in the test section of a flexible-wall wind tunnel as 
illustrated in Figure 1.1a. If the streamlines above and below the model have their 
free air shape, unconfined flow conditions would exist in the test section. These free 
air conditions can be simulated in a wind tunnel if the flow conditions along the 
inner side of the wall are the same as those along the outer side for the fictitious 
external flow field. In principle, this argument is also valid for the ventilated-wall 
tunnels (Figure 1.1b). 
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conventional 


conventional 



adapted 


adapted 


a) Model in flexible-wall wind tunnel 


b) Model in ventilated-wall wind tunnel 


Figure 1.1: A schematic to explain the principle of adaptive- wall concept in 
flexible and ventilated wall wind tunnels. 
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The adaptive-wall wind tunnel provides the means of modifying the boundary 
conditions at the tunnel walls so that the flow can conform to that which would 
exist in free flight. This requires that the flow within the tunnel be made consistent 
at the boundary with a fictitious (computed) flow outside the tunnel that satisfies 
uniform flow conditions at infinity. The tunnel wall conditions must be adapted 
until this consistency is achieved, hence the term ‘adaptive- wall wind tunnel’. 

In practice, the concept of the adaptive-wall wind tunnel requires an iterative 
process in which two flow variables are measured independently in the presence of 
an arbitrary model and flow configuration. At an interface, one flow variable is used 
to evaluate the functional relationships for the unconfined flow. The unconfined flow 
solution is compared with the second flow variable measured at the same interface. 
Once the effect of the wall interference is determined, local wall properties are 
adjusted to achieve the interference-free flow. New flow variables are then measured 
and checked with the unconfined flow solution. If they are not consistent with 
interference-free flow conditions, the procedure is repeated until the results agree 
within the range of the experimental error. 

In general, an adaptive-wall wind tunnel scheme involves three major compo- 
nents: instrumentation for the measurement of flow variables, computing hardware 
and software, and a method of modifying the wall boundary conditions. In ven- 
tilated wind tunnels, adaptive-wall control can be produced by changing the local 
plenum pressure^ 8,9 * or using walls with adjustable porosity^ 10 ). Wall boundary con- 
ditions can be adjusted by changing the local wall contour in flexible-wall wind 
tunnels^ 11 !. 


1.2 Survey of Previous Work 

Since the introduction of the modem adaptive-wall concept, extensive research and 
development has been done on two-dimensional adaptive-wall wind tunnels and 
research on three-dimensional wind tunnels is now in progress. A comprehensive 
bibliography of the studies on adaptive- wall wind tunnels is given in Reference 12. 
The first studies of adapting wind tunnel walls to reduce wind tunnel interference 
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were conducted in the 1940’s* 13 ). Two-dimensional flexible- wall wind tunnels were 
used to solve the blockage problem at high speeds, the shape of the flexible walls 
being determined from the free air solution past the model. This method was limited 
to flows which could be solved analytically. 

The modern adaptive-wall wind tunnel concept was introduced by Ferri and 
Baronti* 14 ) and Sears* 15 ) independently to solve the wall interference problem espe- 
cially in transonic wind tunnels. Ferri and Baronti suggested a method to assess 
wall interference from two independently measured flow variables at the tunnel 
walls. They chose the static pressure and the flow angle as the flow variables. 
These parameters were used in performing free air calculations and the results were 
compared with the measured quantities. In their work, Ferri and Baronti discussed 
the validity of using linearized theory for transonic wall interference assessment and 
demonstrated its applicability by numerical computations. 

In an independent work, Sears developed the concept of a control level removed 
from the walls* 15 ). In his method, two flow variables are measured at a control level. 
One of them is used as a boundary condition to compute the other flow variable in 
unconfined flow external to the chosen control level. The predicted value of the flow 
variable is compared with the measured value to assess the wall interference. An 
iterative procedure is applied to change the wall conditions until the measured and 
computed values agree, implying interference-free flow in the test section. In this 
method, the difficulty of modeling or calibrating the flow near the walls is avoided. 

Application of a numerical simulation of the flow within the wind tunnel test 
section by Erickson and Nenni* 16,17 ), showed that the adaptive- wall wind tunnel 
scheme can be iterated successfully and converges to unconftned flow conditions. In 
this simulation, it was assumed that the flow was incompressible and sufficient wall 
control was available to modify the tunnel wall boundary conditions. 

In the experimental studies performed to show the applicability of Sears’ method 
to two- and three-dimensional flows* 8,18 ), several practical problems related to in- 
strumentation, the nature of the wall modification, and convergence of the control 
logic were discussed in detail. Comparisons with data obtained in larger tunnels 
showed that active wall control reproduced the correct shock wave location and 
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eliminated or minimized the interference effects on aerodynamic characteristics. 

In a theoretical study for the convergence of the adaptive-wall concept, Sears^ 19 ) 
pointed out the major factors influencing the accuracy of adaptive-wall control, 
such as the finite length of the test section and the finite number of control points 
and plenum compartments. Because of these factors, the interference-free rela- 
tionships can be established only approximately. Sears showed that the iterative 
process can still converge to an imperfect approximation to the unconfined flow. 
However, interference can be substantially reduced as compared to that in conven- 
tional, nonadaptive test sections. For most practical purposes, in a real tunnel, 
interference-free flow about the model can be achieved in an adaptive-wall wind 
tunnel with imperfect control. 

For flexible-wall tunnels, compatibility assessment methods based on wall de- 
flection and pressure have been developed and implemented experimentally^ 11,20-23 ). 
To streamline the flexible walls, Goodyear^ 11 ) applied criteria which matches static 
pressures in the test section at a wall with pressures computed for an imaginary in- 
viscid flow field over the outside of the same wall. An iterative procedure is applied 
to change the walls from straight surfaces to streamlines. A very similar adaptive 
test section was designed and installed at NASA Langley for the cryogenic transonic 
tunnel^ 20 ). Several other flexible- wall adaptive tunnels have been designed and put 
into operation for basic and advanced research. Ganzer^ 21 ) reported test results 
from a two-dimensional flexible wall tunnel which was developed with the intention 
of applying adaptive-wall techniques to a test section for three-dimensional model 
tests. 

For subcritical, nonlifting flow, an analytical simulation of an adaptive-wall tun- 
nel by Lo and Kraft/ 24 ) showed that it was possible to find the rate of convergence 
and establish criteria to achieve unconfined flow in a single adjustment regardless 
of the initial conditions. Their analysis was valid for two-dimensional flow fields. 
Dowell^ 25 ) formulated a one-step convergence method which was valid for two- and 
three-dimensional flows. However, his method was subject to the assumption of 
linearity. 
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Adaptive- wall experiments conducted at Arnold Engineering Development Cen- 
ter (AEDC) were aimed at finding out the most suitable wall configuration for an 
adaptive- wall section with particular application to three-dimensional flow^ 10 ). The 
results of the experiments, carried out with different ventilated wall configurations 
with adjustable porosity, showed that wall interference was reduced significantly 
although no completely unconfined flow condition was achieved. 

At Ames Research Center, a two-level one-component compatibility assessment 
was developed by Davis^ 26 h The method is valid for two- and three-dimensional 
flows and was developed to exploit the capabilities of Laser Doppler Velocimetry 
and especially for use with a minicomputer-equipped facility. In this method, only 
one flow variable is measured at two control levels. Values of the flow variable at 
the lower control level provide the boundary condition to solve the unconfined flow 
field at the upper control level. This calculation gives a prediction of the same flow 
variable at the upper control level which is compared with the measured value to de- 
termine the effect of the wall interference. Davis presented numerical simulations of 
two- and three-dimensional flows along with data from a two-dimensional pilot tun- 
nel. The normal velocity component was measured with a one-component LDV 
system and Davis’ two-level compatibility assessment method was implemented 
through the use of influence coefficients to accurately determine the local plenum 
pressure corrections for rapid convergence to interference-free flow^ 9 h 

A new adaptive-wall test section was designed and installed at Ames Research 
Center for a 2x2-foot tunnel on the basis of the experience from the small pilot 
tunnel experiments. It uses a state-of-the art LDV system with a very fast com- 
puter controlled traverse system. Wall control was provided by suction and blow- 
ing through the slotted walls by a number of sliding valves which were operated 
automatically^ 27 ) . 

Development of adaptive-wall test sections for three-dimensional testing has 
been under way for the last decade and is still in the early, demonstration staged 28-31 ). 
At the University of Southampton, a two-dimensional flexible-wall test section was 
used for three-dimensional testing^ 22 ). First results showed that the model blockage 
can be eliminated by streamlining the walls. However, a substantial amount of wall 
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interference remained. It was reported that research was in progress to develop a 
suitable algorithm for the wall interference assessment. 

For the purpose of a very accurate representation of three-dimensional wall 
shape, the test section developed at Deutsche Forschungs und Versuchsanstalt fur 
Luft und Raumfahrt (DFVLR) consists of a cylindrical rubber tube that can be 
deformed by 64 jacks. A one-step method based on the assumption of superposition 
of pressure disturbances, was developed as a calculation procedure to determine the 
adapted wall configuration* 21 ). 

A three-dimensional test section with eight flexible walls was developed at Tech- 
nical University of Berlin* 21,30 ). They form an octagonal cross section and all the 
walls are subject to a two-dimensional deformation. Problems associated with the 
algorithm for the calculation of three-dimensional exterior flow field and the control 
system for wall settings were discussed by Ganzer* 30 ). Test results with full adap- 
tation of the walls demonstrated the feasibility of the adaptive-wall technique for 
three-dimensional testing. 

At Ames Research Center, a two-dimensional test section was modified to per- 
mit cross-stream control for three-dimensional testing* 29 ). The experiments demon- 
strated the feasibility of local wall control in a ventilated three-dimensional test 
section. Although in all the experiments wall interference was reduced, interference 
could not be completely eliminated. It was concluded that the effects of the wall 
adjustments on the flow could not be accurately predicted using influence coeffi- 
cients; in addition, insufficient suction and blowing to produce the required velocity 
changes was provided. 

Application of the adaptive-wall concept to the testing of propulsive lift vehi- 
cles has been theoretically investigated by Sears* 32 ). In a numerical demonstra- 
tion using panel methods, it was shown that high lift configurations can be tested 
in an adaptive- wall wind tunnel of the V/STOL type, correcting for gross errors 
in initial conditions and arriving at approximate unconfined flow conditions. Re- 
cently, limited experimental results were reported from a V/STOL tunnel at Arizona 
University* 33 ). Results showed that the iterative procedure successfully reduced dis- 
crepancies at the control level after about six to eight iterations. 
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1.3 Motivation of the Present Work 

Some of the major problems in the design and operation of ventilated adaptive- 
wall wind tunnels are the complex instrumentation and the excessive total testing 
time. These problems are associated with the type of variable to be measured, the 
data-taking time, and the number of iterations required for convergence. 

Assessment of wall interference and proposed correction schemes for adaptive 
wind tunnels require measurement of two flow variables at control levels at or near 
the walls of the test section. Any two conveniently measured variables can be used 
in the adaptive-wall process. In ventilated test sections, the flow variables can be 
measured using intrusive tools such as probes^ 8 ), and/or nonintrusive tools such as 
Laser Doppler Velocimetry^ 9,35,36 ) (LDV) or side- wall static pressure tappings^ 34- 36 ). 

Because of the wide use of slotted-wall test sections at Ames Research Center, 
adaptive-wall test sections with slotted-walls and local plenum control have been 
developed. Since the side-walls are not porous, the windows can be used for optical 
measurement techniques. At Ames Research Center, LDV is currently used to 
measure velocity components as the flow variables in the compatibility assessment. 
Although it is a nonintrusive and accurate method, collection of velocity data takes 
a long time and the LDV system has to be moved along the control levels to get the 
velocity distribution. It is a complex system and requires special instrumentation. 
Precise optical alignm ent and proper seeding of the flow are necessary to ensure 
accurate data. 

If ventilated adaptive-wall wind tunnels are to be used for production testing 
purposes, adaptive-wall techniques will have to be justified both economically and 
technically^ 37 ). Basically, the adaptive- wall technique is an iterative process and 
total testing time must be reduced for production testing purposes. This provides 
much of the motivation for the present work. 
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1.4 Present Work 

The primary purpose of the present research was to provide solutions to the above- 
mentioned problems by establishing fast schemes with simple instrumentation. In 
order to simplify the instrumentation and to reduce the data taking time, the side- 
wall pressure distribution was proposed as the flow variable^ 34,35 b One of the major 
goals of the present experiments was to validate the use of the side-wall pressure 
distribution in the assessment of wall interference. It is a nonintrusive technique, 
like Laser Doppler Velocimetry, but it is limited to two-dimensional testing. 

As one alternative to reduce the testing time, a one-step convergence scheme 
was developed and implemented theoretically and experimentally. In adaptive-wall 
experiments, each iteration requires the measurement of flow variables and applica- 
tion of pressure corrections to the plenum compartments. The main purpose of the 
one-step scheme is to eliminate intermediate experimental iterations by performing 
them on the computer, using the information from influence coefficient matrices, 
instead of taking the measurements repeatedly. 

Influence coefficients are required to determine the correction to the flow in the 
test section. During actual testing, once the wall interference is determined, the 
amount of suction and/or blowing which would be applied in the plenum chambers 
is determined using influence coefficients. In the present study, influence coeffi- 
cients were obtained in three alternative ways, from side- wall pressure distribution, 
streamwise velocity measurements, and normal velocity measurements. In the itera- 
tive and the one-step methods, influence coefficients were employed for a systematic 
estimation of the pressure corrections. One aspect of the present research was to 
determine the practical limitations and the effectiveness of the influence coefficients 
at different Mach numbers and model incidences. 

One of the purposes of this study was to investigate the relative merits of the var- 
ious convergence schemes with different flow variables. Therefore, the iterative and 
the one-step method were implemented with normal and streamwise velocity com- 
ponents measured using a two-component LDV system and with resultant velocity 
distribution calculated from side-wall pressure measurements. Emphasis was given 
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to the convergence schemes applied with a two-level compatibility assessment^ 26 ) 
which requires the measurement of one flow variable at two control levels unlike the 
one-level method^ 15 ) in which two flow variables have to be measured at one control 
level. 

In this study, linearized equations were used to calculate the unconfined flow at 
the control levels. The linearized compressible potential flow equations are applica- 
ble to supercritical flows past the model, provided that shock waves are not strong 
and do not reach the wind tunnel walls. Therefore, the linearized flow equations 
can be used up to very high subsonic Mach numbers. However, if the flow near 
the walls is nonlinear, the transonic small disturbance equation is the appropriate 
choice to solve the unconfined flow. In this case, the functional relationships can 
not be obtained analytically. In order to investigate the accuracy of the unconfined 
solution for supercritical cases, the transonic small disturbance equation was solved 
using finite difference methods and compared with the unconfined solution obtained 
from the linear equations. 

In Chapters 2 , 3, and 4, theoretical and experimental aspects of the present 
study are explained in detail. Results of the experiments for the evaluation of the 
influence coefficients and the adaptive-wall experiments for the application of the 
various convergence schemes are given in Chapter 5. Results from convergence 
schemes are compared in Chapter 6. Chapter 7 is devoted to the conclusions and 
the recommendations for future work. 


Chapter 2 

Discussion of the Convergence 
Schemes 


2.1 Convergence Schemes 

In adaptive-wall experiments, local wall properties are systematically adjusted to 
achieve unconfined flow by measuring two flow variables and evaluating the func- 
tional relationships for unconfined flow. The iterative process can be entirely ex- 
perimental or it can be carried out computationally by using the one-step conver- 
gence scheme^ 24 ’ 25,35 ). The primary purpose of this study was to apply a two-level 
compatibility assessment with the iterative and one-step convergence schemes and 
determine the relative merits of these two methods. 


2.1.1 Iterative Convergence Scheme 

In the iterative method with the two-level compatibility assessment^ 9,26 ) one flow 
variable is measured at two control levels. A control level is the chosen interface in 
the test section between the tunnel wall and the model, along which the flow vari- 
ables are measured and the compatibility to the free air conditions is assessed. The 
flowchart of the algorithm used for the iterative method is illustrated in Figure 2.1. 
Disturbance velocity components were chosen as the flow variables to explain the 
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steps. After setting the test section Mach number and model angle of attack for 
desired testing conditions, initial settings of the pressures in the plenum compart- 
ments are optional. If they are known approximately prior to the experiments, 
usually from experience, it would possibly help to reduce the number of iterations. 
However, for a different model and Mach number, the logical approach is to start 
with all suction and blowing valves in a closed position. This strategy was chosen 
for the present experiments. plenum 


1 jJLIT jA-L 


CONTROL LEVEL 


MODEL 





Figure 2.1: The flowchart of an algorithm for the iterative method with the 
two-level compatibility assessment 
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The next step in the algorithm is to measure the velocity components at the 
control levels. The measured velocity component at the first level, v mi , is used as a 
boundary condition to evaluate the functional relationships along with the far field 
boundary conditions to determine the unconfined flow solution, v C2 , at the second 
level. The effect of wall interference is determined by comparing the measured 
velocity component with the unconfined solution at the second level. If the measured 
and the calculated velocities fail to agree, then local plenum pressures have to be 
adjusted to eliminate this error and to achieve the interference free flow. Required 
plenum pressure corrections are calculated according to the velocity differences at 
the control points. Since the velocity component measured at the first level is 
not at its interference-free condition initially, the unconfined solution is only an 
estimation, and the pressure corrections will not eliminate the wall interference in 
one step. Instead this operation leads to an iterative process, eventually converging 
to the unconfined flow conditions^ 16,17,19,24 ). After the application of the pressure 
corrections in the plenum compartments, new values of the velocity components are 
measured. If they are not consistent with the interference-free flow conditions, the 
process described above is repeated until the matching of the velocity distribution 
agrees within the range of experimental error. Once this step is reached, measured 
aerodynamic properties and tunnel conditions will be practically free of wind tunnel 
wall interference effects. 


2.1.2 One-Step Convergence Scheme 

The one-step convergence scheme was proposed as an alternative to eliminate the 
wall interference in a single iteration, or to cut down the number of iterations to a 
minimum. The one-step scheme developed and implemented in the present study is 
similar to that proposed by Dowell^ 25 ), using the influence coefficients to determine 
the new distributions of the flow variables at the control levels after adjusting the 
control variables. However, the present scheme is not subject to the assumption 
of linearity and simulates the experimental iterative steps numerically rather than 
trying to compute the required pressure corrections in one calculation. As explained 
in Chapter 5, the elements of the influence coefficient matrices may change during 
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iterations depending on whether suction or blowing is applied. Other one-step 
methods^ 24,25 ) do not have this feature which is necessary in practical applications 
if influence coefficients differ for suction and blowing. The flowchart of the one- 
step convergence scheme is shown in Figure 2.2, where subscripts ‘m’ and ‘p’ stand 
for measured and predicted quantities, respectively. The unconfined flow solution 
is denoted with subscript ‘c’. Although the one-step scheme will be explained for 
the two-level compatibility assessment, it can be easily modified for the one-level 
method. 

After the initial steps, which are exactly the same as the experimental itera- 
tive scheme described in the previous section, measured velocity components on 
the control levels are initialized as the predicted values. Similary, the unconfined 
velocity distribution at the second level is calculated and plenum pressure changes 
required to eliminate the wall interference are found using the influence coefficients. 
In the present method, the pressure changes are assumed to be applied to the plena 
and new velocity distributions at the control levels are predicted using the influence 
coefficients instead of repeating the measurements. These steps are repeated on the 
computer until the predicted velocities, which are the counterparts of the measured 
velocities in the experiments, agree with the calculated values within the prescribed 
error margin. Once this step is reached, the cumulative pressure change for each 
plenum compartment, which is the sum of pressure changes assumed to be applied 
at each step, is calculated. In the actual experiment, these net pressure changes are 
then applied to the active plenum compartments. 

As may be seen from the algorithm, the success of the one-step convergence 
scheme depends on how accurately the velocity distributions are predicted at the 
control levels for given pressure changes. If the number of iterations during the 
experiment were to be reduced to one, the influence coefficients would need to be 
precisely determined. If the influence coefficients are not sufficiently accurate, this 
method may not create a one-step convergence. Even then, the number of iterations 
will be reduced considerably compared to the purely experimental iterative method. 
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In practical applications, the number of control points is usually larger than the 
number of plenum compartments and this constitutes an overdetermined system of 
equations. In the present study, plenum pressure corrections were found in the least 
squares sensed 38 ) for both convergence methods. 


2.2 Influence Coefficients 

The use of influence coefficients, especially in the application of the one-step conver- 
gence method, is the crucial part in the present adaptive- wall experiments. Influence 
coefficients are required to determine the corrections to the flow in the test section. 
For the present ventilated wind tunnel, adaptive-wall control was produced by using 
the localized plenum pressure control. Once the wall interference is determined, the 
amount of suction and/or blowing through the plenum compartments to correct 
this error requires some knowledge of the influence coefficients. Since convergence 
to unconfined flow is an iterative process, the number of iterations depends on the 
accuracy of the influence coefficients as well as the degree of interference. 

An influence coefficient is defined as the change in velocity at a control point for 
a unit change of pressure in the active plenum compartment. If AV| is the velocity 
change produced at control point ‘i’ , and A Pj is the pressure change in the active 
plenum ‘j’, then the slope of AVj vs. A Pj gives the influence coefficient Cij at that 
point. 

The velocity change produced at any control point due to the pressure changes 
in more than one plenum was assumed to be equivalent to the sum of velocity 
changes produced at that point by the same pressure changes in each plenum if 
applied separately^ 9 ) . As a consequence of the assumption of linear superposition, 
it was possible to form influence coefficient matrices from the following equation: 

j=k 

av; = £ c tj AP,, = (1) 

j=l 

where k and 1 are the number of plenum compartments and control points respec- 
tively. Each element of the influence coefficient matrix is defined by C,j = AVJ/APj. 
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Equation 1 can be written explicitly as follows; 



In Equation 1, AVi is the total change in the velocity due to the the pressure 
changes in the active plenum compartments. In the definition of Cij, A Vi is the 
change in the velocity at a control point due to the pressure change in one plenum 
compartment. 

For the construction of the influence coefficient matrix, the velocity change at 
a control point was obtained using the following algorithm. The pressure in each 
plenum compartment was varied systematically and the desired velocity component 
was measured at each point on the control levels 1 and 2 which are illustrated in 
Figure 2.3. For each plenum compartment, initial velocity distributions were taken 
with all the control valves in closed positions. In the subsequent runs, suction 
or blowing was applied in the active plenum compartment at five different valve 
settings while the control valves for the rest of the plenum compartments were kept 
closed. This procedure was repeated for each plenum compartment. 

The velocity change at a control point due to the pressure change in a plenum 
compartment was obtained using the following equation: 

AVi = ( U , - U Tef ) (n) - (Ui - U ref \ o) (3) 

where i is the i th control point and U re f is the freestream velocity used as a 
reference value. Ui is the velocity calculated from side-wall pressure distribution or 
the streamwise velocity component measured with the LDV system. In Equation 3, 
‘n’ represents the case when suction or blowing is applied in the active plenum 
compartment, and ‘o’ refers to the case when all valves are in the closed position. 

For the normal velocity component, the reference velocity is assumed to be 
negligible upstream, since the normal disturbance velocity component goes to zero 
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far away from the model. Thus, the change in the normal velocity at a control point 
can be written as 


AV; = (u,-)( n) - (v, •)(<>) (4) 

The pressure change in a plenum compartment is simply taken as 

A Pi = (P;) w - (Pj) W (5) 

where j is the j th plenum compartment. 

Since the plenum compartments were arranged symmetrically with respect to the 
centerline of the test section, they could be grouped in pairs. In this work, unless 
otherwise stated, plenum compartment j refers to both top and bottom plenum 
compartments located symmetrically at the j th location (See Figure 2.3). During the 
experiments carried out to determine the influence coefficients, suction or blowing 
was applied symmetrically through top and bottom plenum compartment j. 

A series of experiments was carried out to evaluate the influence coefficients from 
the side- wall pressure distribution, streamwise velocity distribution, and normal ve- 
locity distribution^ 34,35 ). For each plenum compartment and control point, influence 
coefficients were evaluated for suction and blowing at different Mach numbers with 
and without the model. In the case of the side- wall pressure distribution, influence 
coefficient matrices were formed by the data taken without the model, while those 
for the streamwise and the normal velocity components were constructed from the 
data taken in the presence of a nonlifting model. The effects of model, Mach num- 
ber, location of control points, suction and blowing on the influence coefficients 
were evaluated. Linearity and the validity of superposition of the influence coeffi- 
cients were investigated. Experimental results and the evaluation of the influence 
coefficients are discussed in detail in Chapter 5. 
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Figure 2.3: Test section arrangement for the adaptive-wall tunnel experiments. 



Chapter 3 

Unconfined Flow Analysis 


Calculation of the unconfined flow is the principle theoretical aspect of the adaptive- 
wall concept. A suitable functional relationship is evaluated to find the unconfined 
flow field external to the chosen control level. 

In the present experiments, control levels were far enough away from the model 
and top and bottom wall boundary layers so that the regions of complicated flow 
patterns and viscous effects such as separated flow regions, shock-boundary layer 
interaction, were avoided. The undisturbed freestream velocity was always subsonic 
and the disturbances at the control levels were assumed to be sufficiently small that 
the small disturbance equations would be a valid approximation. Throughout this 
study, linearized compressible potential flow equations were used to solve the uncon- 
fined flow field. To employ the more complex nonlinear transonic small disturbance 
equations on line, a powerful dedicated computer would have been required. 


3.1 Linear Solution 

The governing partial differential equation for two-dimensional linear subsonic flow 
is the Prandtl-Glauert equation which can be expressed as: 

0 2 <t>xx + 4>yy = 0 ( 6 ) 
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where 

P = y/i- Moo 2 


Equation 6 was solved in two different ways in accordance with the type of 
compatibility assessment used in the adaptive- wall experiments. 

For the one-level two-component assessment method, the solution of the un- 
confined flow field is related to the solution of Equation 6 for a half space using 
Cauchy’s Integral Equation^ 39 ) (Figure 3.1). The unconfined flow solution for the 
normal velocity component at y = ±h, for the boundary conditions given in terms 
of the streamwise disturbance velocity, along with the far-field condition that the 
streamwise disturbance velocity vanishes, is given as^ 7,15 h 


v(x, ±h) = ±— 4) 

7T J-oo 


P r + °° u m (Z,±h)d£ 


( 7 ) 


y . *n 


z-plane 



Figure 3.1: Computational domain used for the calculation of unconfined flow 
using Cauchy’s Integral Equation. 
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If the boundary conditions are given in terms of the normal disturbance velocity, 
then the unconfined flow solution for the streamwise disturbance velocity component 
is: 


_ _ 1 r + °° VmU,±h)dt 


u{x ' ±h) = *~pK i 


x 


( 8 ) 


Equations 7 and 8 must be evaluated in the sense of the Cauchy principal value. 
For the present applications, the singularity in each of these equations was removed 
using a method suggested in Reference 40 and the integral was evaluated by using 
spline fit^ 41 ) values of the measured velocities, v m or u m . The derivation for the 
removal of the singularity is given in Appendix A. 

For the two-level one-component assessment method, a solution to the Prandtl- 
Glauert equation for a half plane using the normal velocity component was sug- 
gested by Davis^ 26 ). By differentiating Equation 6 with respect to x, the following 
differential equation for the streamwise velocity component can be obtained: 


P 2 u xx + u yy = 0 (9) 

The following boundary conditions are imposed to solve the unconfined flow field 
in the region shown in Figure 3.2: 

u = u m (x) at y = hi and — oo < x < +oo 

u — > 0 as r = y/x 2 -f y 2 — ► oo 

Solution to this boundary value problem can be found using Poisson’s Integral 
Formula for a half plane^ 39 ) along with the Prandtl-Glauert correction. For the 
streamwise disturbance velocity component, it can be written as: 


u(x,y) = 


P(y~hi) f +0 ° 


7 r 


I- 


(*- 0 2 + /? 2 (y-M 2 

and at the upper control level, y = h 2 , it is simply as follows: 


, L ^ P h r +°° «„»(£,/»!)<*£ 

= (x _ {)2 + J, v 


+oo 


( 10 ) 


(ll) 


where h = h 2 — h\. The same solution is also valid for the lower region. Since the 
boundary conditions are imposed at the control levels independently in each region, 


CHAPTER 3. UNCONFINED FLOW ANALYSIS 


24 


region, the flow field is decoupled with the requirement that the measurements must 
be made along the boundaries far enough away from the model. This is necessary 
to insure that the contributions from both regions are small. It is noted that the 
integral in Equations 10 and 11 is nonsingular for the two-level compatibility method 
compared to the one-level assessment method. It can be evaluated using a simple 
numerical integration scheme. y 



Bottom wall 

Figure 3.2: Computational domain used for the calculation of unconfined flow 
for two-level assessment method. 


3.2 Nonlinear Solution 

It is appropriate to use the transonic small disturbance equation when the flow is 
supercritical at the control surfaces. However, during the experiments, the com- 
puting power of the dedicated computer was not sufficient to employ it on line. 
Later, the solutions obtained from the linear equations were compared with the 
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i solutions from the nonlinear equation to determine the accuracy of the method for 

the present Mach number range. 

The governing equation for the small disturbances in the transonic flow field is 
! given as( 4 h 

+ = 0 ( 12 ) 

U oo 

where <j> x = u and <f> y — v are the streamwise and the normal disturbance velocity 
components respectively. 

Boundary conditions for the streamwise velocity component are imposed as fol- 
lows: 

(f> x = u — > 0 as x and y — > oo 

<f> x = u m (x) at y = ±hi and for — oo < x < -f oo 

Equation 12 was solved numerically using the Murman-Cole method with Gauss- 
Seidel line relaxation in the computational plane shown in Figures 3.3a and b. The 
derivation is given in Appendix B. 

y 



Figure 3.3a: Physical domain with boundary conditions for the transonic small 
disturbance equation. 
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112 x 36 computational mesh points 

Figure 3.3b: Computational mesh for the calculation of unconfined flow using 
transonic small disturbance equation. 

3.3 Discussion 

As previously stated, one of the purposes of this study was to use the resultant 
velocity determined from the side-wall pressure distribution as the flow variable. 
Since the flow angle was not measured in the present experiments, the calculated 
resultant velocity distribution could not be resolved into its streamwise and nor- 
mal components. The velocity measurements taken with the LDV system showed 
that the magnitude of the resultant velocity was comparable to the streamwise 
component. For the resultant velocity distribution from side-wall pressure mea- 
surements, unconfined flow calculations had to be carried out approximately by 
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assuming that the streamwise velocity component is large compared to the normal 
disturbance velocity component. In order to find the error introduced with this 
assumption and the feasibility of this approach in the adaptive-wall concept, the 
iterative steps of the schemes with the streamwise velocity distribution were recon- 
structed computationally implementing the resultant velocity distribution in the 
unconfined flow calculations. In this case, resultant velocity distributions were cal- 
culated from streamwise and normal velocity components which were measured with 
the LDV system. Results and discussion of this assumption are given in Chapter 5. 

Since the flow is symmetrical with respect to the centerline of the test section 
at zero angle of attack, only the upper half of the computational plane was used 
to solve the unconfined flow field. In the case of a lifting airfoil, unconfined flow 
calculations were performed separately for the upper and the lower computational 
planes. It was assumed that the upper and lower flow fields were decoupled for 
small model incidences. 


Chapter 4 

Experimental Apparatus and 
Data Acquisition 


4.1 Experimental Setup 

4.1.1 Wind Tunnel 

Experiments were conducted in the ventilated adaptive-wall wind tunnel at Ames 
Research Center. It is a continuous flow, indraft wind tunnel equipped with a flex- 
ible downstream throat to control the test section Mach number, and an auxiliary 
system to provide suction and blowing of air through the plenum chambers. The 
adaptive-wall test section with slotted walls is of rectangular cross section, 25 cm. 
x 11 cm. and is 61 cm. long. The tunnel is operated at atmospheric pressure. The 
upper and lower walls each have 10 longitudinal slots with a total open area ratio 
of 12-percent. The test section has separate top and bottom plenums and each 
was divided into 10 plenum compartments by transverse inserts. The transverse 
plates can be removed and inserted easily, so that number of plenum chambers can 
be varied in increments of 2.54 cm. A schematic of the wind tunnel along with a 
typical connection of a plenum compartment to the auxiliary system is illustrated 
in Figure 4.1. 

The tunnel is operated with the throat choked to assure a steady flow in the test 
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section and to minimize the noise propagated from the control gate valve. Throat 
height is changed by using an electric motor whose speed is controlled by a drive 
unit. A digital voltmeter is used to monitor the throat height. 



Figure 4.1: Schematic of the 25 cm x 11 cm NASA Ames indraft wind tunnel. 
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Each plenum compartment was connected to an auxiliary air system which pro- 
vided suction and/or blowing when required (Figure 4.1). Blowing in a plenum 
compartment was obtained by opening the corresponding pressure valve to atmo- 
spheric pressure. An ejector which provided 16 inches of Hg. vacuum pressure 
(gage) at the maximum setting, was used for suction. It was connected to a vac- 
uum manifold to which all the plenum compartments were connected. During the 
experiments the vacuum and pressure control valves were operated manually. 

4.1.2 Model 

For the experiments, a 7.62-cm.-chord NACA 0012 airfoil without boundary layer 
trips was used. It had a 25 cm.-span and was instrumented with 16 pressure tap- 
pings on the upper surface and 9 tappings the lower surface. The blockage ratio 
of the model in the test section was 8.3 % and it was expected to create signifi- 
cant wall interference effects at high speeds. The aspect ratio of the model, tunnel 
width/model chord, was 3.28 and tunnel height/model chord ratio was 1.44. This 
airfoil section was chosen because of the large body of experimental data and numer- 
ical solutions that exist for its aerodynamic characteristics, especially at transonic 
speeds. It is used for helicopter rotors and has been selected as one of the standards 
by AGARD to be used as an experimental data base for the validation of advanced 
Computational Fluid Dynamics codes. Figure 4.2 shows the airfoil section with 
orifice locations. 


4.2 Instrumentation 

Figure 4.3 shows the block diagram for the instrumentation used in the adaptive- 
wall experiments. A pitot tube and thermocouple unit were mounted upstream of 
the test section to measure the total pressure and ambient temperature. 
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Figure 4.2: Pressure orifice locations of NACA 0012 airfoil model. 


4.2.1 Instrumentation for the Pressure Measurements 

Plenum, model and control level pressure distributions were measured using Scani- 
valves. Calibration pressure was supplied by a vacuum pump and measured with 
a separate unit. A Paroscientific pressure computer was used to measure ambient 
pressure. Throat, test section and plenum pressures were also monitored on the 
manometer board during the experiments. Side-wall static pressure distributions 
at the control levels were measured with a number of pressure tappings located on 







CHAPTER 4. EXPERIMENTAL APPARATUS AND DATA ACQUISITION 32 


one of the side-walls. A six-unit Scanivalve system was used to measure the stag- 
nation pressure, top and bottom plenum pressures, top wall pressure distribution, 
downstream throat pressure and side wall pressure distribution. 



Figure 4.3: Adaptive- wall wind tunnel instrumentation. 
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Each Scanivalve unit had 24 ports which were connected to the pressure tap- 
pings. The side-walls were made of plexiglass and one of them was instrumented 
with seven rows of pressure orifices along the test section with a total of 141 pres- 
sure tappings. They were located symmetrically above and below the centerline. In 
order to best resolve the flow field, the pressure tappings were closely spaced near 
the region where the model was located (Figure 4.4). 

Calibration of the Scanivalves 

The transducers in the Scanivalve system were of the differential type with a range 
of ±5 psig. In order to calibrate them, a calibration pressure generated by an electric 
vacuum pump was applied in one port of each Scanivalve unit. Another port was 
tapped to the atmosphere. The atmospheric and calibration pressures were each 
measured by alternate methods. These two known pressures and the values given 
for them by each transducer allowed the Scanivalves to be calibrated. 
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Figure 4.4: Adaptive- wall wind tunnel test section arrangement for the side-wall 
pressure measurements. 
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4.2.2 Laser Doppler Velocimeter 

Velocity measurements at the control points were taken with a two-component Laser 
Doppler Velocimeter system. A schematic of the LDV system is given in Figure 4.5. 

A 5W Argon-Ion laser (Spectra-Physics Model 64) was used in the forward- 
scatter LDV system. The laser is configured with an end cavity multiline reflector 
which allows it to operate in the multiline mode. In this operating mode, the 
resultant beam is comprised of all wavelengths which are inherent to the Argon- 
Ion uni t.. An external collimater is not required since the resultant beam exits the 
laser in a fairly well collimated state. A broadband polarization rotator mounted 
at the exit of the laser is used to rotate the polarization vector from the vertical 
to the horizontal plane before the color separator unit which comprises of a pair 
of Brewster angle prisms. The polarization of the beam should be in the plane 
formed by the beam and the normal to the prism surface to obtain the maximum 
transmittance* 42 ). Two higher power color lines, green (514.5 nm) and blue (488 
nm), are guided by means of turning mirrors to the Bragg cells. Before the Bragg 
cell, the green beam is passed through a half-wave plate to rotate the plane of 
polarization 90° so that all the optical components for the green line will lay within 
the same plane to avoid the unwanted rotations of the polarization vector. Each 
color line is split into two equal-intensity, slightly diverging beams while they are 
passed through the Bragg cells. A Bragg cell is an acousto-optical device which 
induces a 40 MHz modulation in one of the beams in each pair. For each color, 
the two diverging beams are then made parallel by passage through an alignment 
prism. The prism should lay in the same plane with the Bragg cell to insure that 
the parallel beams will not be skewed. Small corrections are made possible using a 
polarizing beam-splitter cube while filtering out the unwanted polarization. Turning 
mirrors are used to guide the resulting beams through the translation stages. By 
means of a lens, the beams are crossed at the desired point in the test section. Prior 
to the lens, all four beams must be parallel and symmetrical to the axis of the lens 
to focus to a common point. 

The receiving optics is an off-axis, forward-scatter configuration which gives 
flexibility to the system. It can be adjusted for different requirements; if maximum 
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signal strength is required, angular deviation from the axis of symmetry should be 
kept to a minimum. Angular deviation can be kept as large as possible, depending 
on the physical limitations, if focusing on the smallest part of the probe volume is 
of interest. The first lens in the receiving optics collects the scattered light and the 
second lens makes the collected beam parallel. A parallel beam is required by the 
polarization beam splitter and the line filters for maximum efficiency. While the 
green lines pass straight through the beam splitter, blue lines are deflected upward. 
The line filters eliminate any cross-talk caused by impure polarization. Each line 
filter is followed by a lens which focuses the light onto a pinhole aperture of the 
corresponding photomultiplier tube (PMT). The PMT’s convert the light energy 
to electrical signals which are passed on to the signal processors. Preamplifiers are 
located in the PMT housing to increase the signal amplitude before driving relatively 
long cables. The preamplifier is kept close to the photodetector in order to increase 
the signal level before additional electronic noise affects the signals. An additional 
amplifier is also used to increase the signal strength prior to the high pass filter. A 
sine- wave generator provides the mixing frequency before the passive diode mixer, 
after which the output is fed to the signal processors. The signal processors used 
in the experiments were of the counter type (Microdyne 2000 series) which have 
good resolution and noise handling characteristics. These units employ circuitry 
with bandwidths in excess of 1 GHz which has been applied under single particle, 
multiparticle and multivelocity conditions. In addition to the basic 5/8 (or 10/16) 
signal validation criteria and oversize particle rejection, a sophisticated multilevel 
sequence validation has been added. This allows the experimenter to select the 
criteria which will yield best data for a given signal-to-noise ratio, signal level, 
seeding level, doppler frequency and fluctuation, indicative of the test situation. 

Calibration of the Laser Doppler Velocimeter System 

At the probe volume, the fringe spacing, 6, is determined using the following 
formula^ 42 ) ; 
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2 sin( |) 

where A is the wavelength and 6 is the angle of the intersecting beams. 

Knowing the fringe spacing and the time required for a given number of crossings, 
particle velocity can be determined. In practice, this method may not give accurate 
results because of small misalignments of the optics. It is best to calibrate the 
system by measuring the scattered light from an object moving at a known velocity. 
A rotating disk is usually the most convenient tool for this purpose. 

TWO-COMPONENT LASER DOPPLER VELOCIMETER SYSTEM 



PHOTOMULTIPLIER 



Figure 4.5: Schematic of the Laser Doppler Velocimeter system. 
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Seeding Requirements 

Seeding of the flow is an inconvenient necessity of the LDV technique. The LDV 
actually measures the velocity of the small particles which are assumed to move 
with the flow. Size of the particles and the concentration distribution are important 
factors which may adversely affect the results. For most gas flows, particle sizes 
on the order of 1 fim are often accepted as adequate^. Smaller particles will 
contribute to the background noise although they may not be detectable by the 
instrument while larger particles may not adequately track the flow. 

In the present study in an attempt to reduce the data taking time, the flow 
was seeded with 0.5-micron-diameter polystyrene spheres diluted with denatured 
alcohol . Analysis and the experimental results^ 44 ) have shown that these particles 
respond to the motion of aerodynamic test flows with excellent accuracy. This 
seeding technique also avoids deposition of the seed material on the test section 
windows as is the case for oil-based seed. 


4.3 Data Acquisition 

The number of control points was variable in the application of the convergence 
schemes with the maximum number being 15 control points at each control level. 
A minicomputer dedicated for the adaptive-wall wind tunnel experiments was used 
for the data acquisition and real-time data analysis. Pressure, velocity and tem- 
perature measuring devices were interfaced with the analog-to-digital converter of 
the computer. A traversing system, whose driving motors were controlled by a pro- 
grammable control unit interfaced with the computer, was used to move the probe 
volume to the desired location for the LDV measurements. 

4.3.1 Velocity Calculations from Pressure Measurements 

Assuming isentropic flow, the Mach number, M can be calculated from the stagna- 
tion pressure, P t , and the static pressure, P, measured at the control point: 
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M = 

Then, the velocity is given by 




(14) 


V = Ma (15) 

where 

a = y/yRT (16) 

where 7 = 1.4 and T is obtained from the total temperature, T t , using the 
following formula: 


T = 


T t 



(17) 


4.3.2 Velocity Calculations from LDV Measurements 

The Macrodyne signal processor used in the experiments utilizes a 1 GHz counter. 
The number of counts for the passage of 8 cycles is sent to the digital I/O in the form 
of a 10-bit mantissa and 4-bit exponent. The signal frequency, f s . is determined 
from the following equation ( 43 h 


1 2 5 "’ exponent 1000 

t s mantissa 

where f s is obtained in MHz and t, is the period of the signal when 8 cycles are 
measured. Exponent and mantissa are set on the signal processor for the desired 
signal conditions. 

Since the signal frequency, / a , is a linear combination of the Doppler frequency, 
fo, the shift frequency, /shift-, and the mixing frequency, / mtx , the Doppler frequency 
can be obtained from: 



f D — fs f shift “b fmix 


(19) 



! 
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Once the Doppler frequency of a particle is known, its velocity can be determined 
from: 

V = Co/ D (20) 

where c 0 is the constant which is determined based on the fringe spacing or 
direct calibration. 

4.3.3 Reynolds Number and Model Pressure Coefficient 

Test section Reynold Number based on the airfoil chord length is calculated from: 

Re = — (21) 

v 

where and c are the freestream velocity and the model chord length respec- 
tively. v denotes the kinematic viscosity. 

Model pressure coefficients are calculated from the measured pressure distribu- 
tion on the model surface from: 


P~ Poo 


Cp Wvi 

or by using compressible flow equations, Equation 22 is arranged as: 


( 22 ) 


Cp = 


(r)(i + 


iMl IV'- ' 2 

where p/p t and M <» are obtained from measurements. 


(23) 


4.4 Error Analysis 

4.4.1 Error Analysis in the Measurements 

Pressure Measurements 

An average of 200 samples in a period of 5 microseconds was taken as the mean 
value in the pressure measurements. Precision of the A/D converter is given as ±1 
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count in the full range of 16383 counts at maximum. Repeatibility of the pressure 
measurements were accurate up to the last digit in counts which corresponds to a 
relative error in a range from 0.1 % to 0.05 %. 

Velocity Measurements 

For velocity measurements, 1000 samples were taken by the LDV system in each 
velocity direction (x and y) at every measurement location. The coincidence of the 
data from the two channels was confirmed within the controlling software. In the 
data reduction, statistical mean and standard deviation based on 1000 samples were 
computed first. A histogram was formed plotting the frequency versus the number 
of samples at each frequency to determine if the data is conforming to the normal 
distribution as a real time performance check. Standard deviation was of the order 
of 1.0 m/s at lower Mach numbers (M ~ 0.5). At higher Mach numbers (M ~ 0.7), 
it increased slightly, but it stayed less than 2.0 m/s. The average velocity at each 
location was calculated based on the number of samples that fell into a region 
bounded by a value of the two-thirds of the maximum value on the histogram. This 
procedure was used to improve the data quality by removing any bias error. In 
the experiments, the number of samples used in the velocity calculations with this 
method varied between 700 and 950. Mean velocity based on 1000 samples usually 
agreed with the velocity calculated from the new procedure of the order of ~ 0.5 m/ s 
at lower Mach numbers. At higher Mach numbers, this value increased, possibly 
because of condensation effects. It should be noted that the standard deviation 
values given here were based on 1000 samples. Therefore, it was concluded that the 
new procedure gave better estimates. Repeatibility of the velocity measurements 
varied between ±0.3 m/s and ±0.6 m/s on average. 

4.4.2 Error Analysis for the Convergence 

In the present experiments, the root-mean-square (rms) error criteria was used 
as a measure for the convergence to the interference-free flow conditions in the 
applications of the adaptive-wall schemes. The rms error in terms of streamwise 
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velocity was defined as; 


rms = 


l 1 

^ J — U c,i ) 2 


(24) 


where u m ^ = U m ^—U re j and u Ct , is the calculated disturbance velocity component 
at i th control point. ‘1’ denotes the number of control points. The rms error was 
similarly defined for the normal velocity and the resultant velocity from the side- wall 
pressure. 


4.4.3 Effect of the Side- Wall Boundary Layers 

Side- wall pressure measurements were used to compute the corresponding resultant 
velocity distributions at the control levels. The presence of the side-wall boundary 
layers was found to have negligible effect on the velocity calculations during the 
test section calibration. The maximum variation in the Mach number measured 
upstream and downstream of the test section was ±0.002 at M — 0.84 and ±0.001 or 
less at lower Mach numbers. These values translate to a variation of ±0.5 m/s or less 
in the velocity calculations in the Mach number range of the present experiments. 
At M = 0.7, velocity measurements taken by the LDV system at the far-field 
control points upstream and downstream of the model showed a deviation of ~ 1.0 
m/s which was found to be in agreement with the side- wall pressure measurements. 

Since the present slotted adaptive-wall tunnel has a different test section ar- 
rangement (plenum compartments, suction and blowing), previous studies on the 
boundary layer growth and its effects on the wind tunnel walls in conventional 
slotted-wall wind tunnels, were not used. However, as explained in the previous 
paragraph, calibration data did not show any unusual flow behavior in the test 
section. 
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4.5 Calibration of the Wind Tunnel 

The relationship between the test section Mach number and the flexible throat 
height was established by calibration. The test section Mach number was set ap- 
proximately using this calibration, however, precise Mach number adjustments were 
made in accordance with the pressure measurements. 

The calibration of the test section was carried out at different Mach numbers 
without the model and with all the plenum valves closed. Centerline Mach number 
distributions along the test section were obtained from pressure measurements on 
the side-wall. In Figure 4.6, the distributions are shown for M = 0.3 to M = 
1.0. Flow develops smoothly in the convergent nozzle at the entrance of the test 
section. Downstream, the cross-sectional area of the test section grows slightly to 
accommodate the boundary layer growth. However, at M = 0.98, this viscous effect 
dominates and the flow chokes before the downstream throat. 



x/c 

Figure 4.6: Centerline Mach number distribution. 
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The test section Mach number distribution in the axial direction was determined 
between X/C = —3.0 and X/C = 3.0 . Up to M = 0.85 , the standard deviation 
from the mean is less than or equal to ±0.003 and it is a reasonably accurate refer- 
ence value. Figure 4.7 shows the test section Mach number distribution along the 
axial direction. Variation of the Mach number distribution in the vertical direction 
was also investigated at various locations upstream and downstream of the model 
location. Typical distributions obtained at two different x-locations are shown in 
Figure 4.8. 



Figure 4.7: Test section Mach number distribution in axial direction. 
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Y/C 


Figure 4.8: Test section Mach number distribution in vertical direction. 

The streamwise gradients in the Mach number distribution can easily be elimi- 
nated by proper blowing or suction. As shown, a reasonably good distribution was 
obtained up to M = 1.0 in the empty test section. At freestream Mach numbers 
above 0.98, the shock moved forward into the test section and operation of the 
downstream throat became ineffective. Further tests done with the model in place 
showed that a maximum Mach number of 0.75 could be reached with the throat 
choked. Higher throat openings did not produce any change once this maximum 
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value was reached. 

If the entrance length is not long enough, it is always hard to determine the 
test section Mach number in the presence of the model. In order to overcome this 
difficulty, the pressure in the plenum compartment 1 was chosen as the reference to 
determine the freestream properties. The plenum calibration procedure is consid- 
ered as a standard method of defining the test section Mach number by relating the 
plenum pressure to an average test section centerline Mach number. The relation- 
ship between the test section Mach number and the Mach number corresponding to 
the state in the first plenum was determined. Because of the size and the location 
of this compartment, it was thought to be a reasonable choice to determine the 
freestream velocity and Mach number for the adaptive-wall experiments. In the 
presence of the model, velocity measurements taken upstream of the model were 
compared with this reference for calibration. In Figure 4.9, freestream velocity and 
Mach number calibration curves are given for the reference plenum compartments. 



Figure 4.9: Calibration curve for freestream velocity. 


Chapter 5 
Experiments 


Adaptive-wall experiments were carried out in two phases. First, influence coeffi- 
cients were evaluated systematically for different flow variables and Mach numbers. 
In the present study, evaluation of the influence coefficients for application to the 
convergence schemes is very crucial to accurately determine the plenum pressure 
corrections. 

In the application of the convergence schemes, emphasis was given to those 
using the streamwise velocity component and the resultant velocity from side- wall 
pressure measurements. Experimental studies with the normal disturbance velocity 
component were completed in a similar test section^. 

In these experiments, the minimum spacing for the plenum compartments was 
over the model, a region where the velocity gradients are high. Eight out of ten 
plenum compartments were used. The first plenum compartment was chosen as the 
reference plenum for the freestream values and the last plenum compartment was 
shown to be ineffective for control of the flow. 

Velocity measurements were taken at the centerline of the test section along 
the streamwise direction. When the LDV system was used, the side-wall with the 
static pressure tappings was replaced by a clear side-wall. The number of control 
points and their locations were the same as those used in the side-wall pressure 
measurements. Measurement of the velocity components at the control points by 
the LDV system consumed most of the available time. 
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5.1 Evaluation of the Influence Coefficients 


Experiments for the evaluation of the influence coefficients were carried out in two 
parts. In the first part, the static pressure distributions at the control levels were 
measured by a number of pressure tappings located on one of the side-walls. The 
maxim um number of control points at each level was 18 although the influence 
coefficients were evaluated at 15 points. In the second part, using the LDV system, 
streamwise and normal velocity components were measured at 8 points at each 
control level to determine the influence coefficients. Because of time limitations, the 
number of control points was kept at the minimum. Velocities at the intermediate 
points were calculated by cubic spline interpolation. To determine the effects of 
model and Mach number, measurements were taken with and without the model 
and at different Mach numbers. Model angle of attack was kept nominally at zero 
degree. The test conditions are summarized in Table 1. 


Side-wall pressure measurements without model 

M 

0.5 

0.7 

0.8 

Test 

Condition 

Plenums from 1 to 10 
Plenums 4 and 6 comb 
Top and bottom pi. 6 

Plenums 4 and 6, 
Plenums 4 and 6 
combined 

Plenums 4 and 6, 
Plenums 4 and 6 
combined 

Side-wal 

pressure measurements with model, a=0 deg 

M 

0.5 

0.685 

0.75 

Test 

Condition 

Plenums 4 and 6, 
Plenums 4 and 6 
combined 

Plenums 4 and 6, 
Plenums 4 and 6 comb 
Top and bottom pi. 6 

Plenums 4 and 6, 
Plenums 4 and 6 
combined 

LDV measurements without model 

M 

0.5 

0.7 


Test 

Condition 

Plenums from ,2 to 9; 
Plenums 4 and 6 
combined 

Plenums 4,6,8; 
Plenums 4 and 6 
combined 


LDV measuremen 

ts with model, a=0 deg 

M 

0.5 

0.7 


Test 

Condition 

Plenums from 2 to 9; 
Plenums 4 and 6 
combined 

Plenums 4,6,8; 
Plenums 4 and 6 
combined 



Table 1: Test conditions in the experiments for the evaluation of the influence 
coefficients. 













CHAPTER 5. EXPERIMENTS 


48 



5.1.1 Effect of Suction and Blowing 

Resultant velocity distributions at the control levels were calculated from side-wall 
pressure distributions using compressible flow relations (Chapter 4). Streamwise 
and normal velocity components were measured using the LDV system described in 
the previous chapter. Experiments showed that the disturbance velocity distribu- 
tion at the control levels could be altered significantly by varying the local plenum 
• 

pressures. 

In Figure 5.1a , the change in the normal velocity distribution at level 2 is shown 
when different suction and blowing pressures were applied in plenum 4 at M = 0.5 
without the model. In the presence of the nonlifting model, the normal velocity 
distribution at level 2 is altered by suction and blowing in plenum 7 as presented 
in Figure 5.1b. For the normal velocity component, the maximum changes were 
produced at control points immediately below the active plenum compartment. 
Suction and blowing produce velocity components toward and away from the wall, 
respectively. 

PLENUM NO. 4 M = 0.5 LEVEL 2 WITHOUT MODEL 



-3-2-10 1 2 3-3-2-101 23 

X/C X/C 


Figure 5.1a: Variation of the normal disturbance velocity distribution with suc- 
tion and blowing. 




Figure 5.1b: Variation of the normal disturbance velocity distribution with suc- 
tion and blowing. 

Compared to the normal velocity component, a completely different velocity 
distribution was observed for resultant and streamwise velocities. As illustrated 
in Figure 5.2a, the effect of plenum 4 on the streamwise disturbance velocity at 
points upstream of that plenum is negligible compared to downstream points. In 
Figure 5.2b, the disturbance in resultant velocity distribution at level 2 shows the 
same behavior at a different Mach number when suction and blowing are applied 
in plenum 4. The same trend is observed for plenum 6 in the presence of a model 
and at a different Mach number as shown in Figure 5.2c. The effect of suction 
is to produce a negative disturbance velocity while blowing produces a positive 
disturbance. Figure 5.3 illustrates the behavior of the influence coefficients for the 
resultant velocity at various control points along levels 1 and 2 for plenum 4 at 
M = 0.5. 
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Figure 5.2: Variation of the streamwise and resultant velocity distributions with 
suction and blowing. 
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PLENUM No. 4 M = 0.50 WITHOUT MODEL 

□ LEVEL 1 SUCTION * LEVEL 1 BLOWING 
a LEVEL 2 SUCTION v LEVEL 2 BLOWING 











-.04 -.02 0 .02 -.04 -.02 0 .02 -.04 -.02 0 .02 

AP/P t AP/P t AP/P t 


Figure 5.3: Influence coefficients at various control points for resultant velocity. 
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5.1.2 Effect of Mach Number 

In general, the dependence of the influence coefficients upon the Mach number and 
presence of the model is as expected. The effects of Mach number and model on the 
influence coefficients are crucial for practical applications. The degree of dependence 
and the limitations should be carefully explored. 

The effect of Mach number on an influence coefficient at control point 12 when 
suction and blowing is applied into plenum 4, is illustrated in Figure 5.4a. Veloc- 
ities were calculated from side-wall pressure distributions which were measured in 
the presence of the nonlifting model at different Mach numbers. For blowing, the 
influence coefficients change very little with Mach number. However, differences 
are observed in the suction case. The effect of plenum 4 is negligible at upstream 
control points for the Mach number range of interest. 



Figure 5.4a: Effect of Mach number on the influence coefficients evaluated from 
side-wall pressure distribution. 
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Figure 5.4b: Effect of Mach number on the influence coefficients evaluated from 
side- wall pressure distribution. 

Figure 5.4b illustrates the effect of Mach number at control point 12 for plenum 
6 with the model present. At M = 0.75, the influence coefficient shows a clear 
departure from the others in the blowing case. The location of plenum 6, which 
is directly above the model, suggests that blowing may contribute to a shock wave 
interaction which causes this discrepancy. In Figure 5.4a, this anomaly was not 
observed for plenum 4, which is located upstream of the model. 
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Figure 5.5a shows the influence coefficients with the model present evaluated 
from the normal velocity distribution at control point 5 which is right below plenum 4. 
Figure 5.5b shows influence coefficients for the streamwise component without the 
model created by about the same blowing and suction. The effect of Mach number 
is negligible for both velocity components. 

The differences at some points on the same curve are of the order of 1 m/s which 
is about the accuracy of the LDV measurements. Velocity measurements with LDV 
are very sensitive to alignment errors and non- uniformities in the plexiglass side- 
walls. 



Figure 5.5a: Effect of Mach number on the influence coefficients evaluated from 
normal velocity distribution. 
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Figure 5.5b: Effect of Mach number on the influence coefficients evaluated from 
streamwise velocity distribution. 

5.1.3 Effect of Model 

Another important aspect for the influence coefficients which should be explored 
carefully, is the effect of. the model. If they are dependent on the type of model 
to be tested, the tedious and costly determination of the influence coefficients has 
to repeated for different models. The same argument is also valid for the Mach 
number effects. 

Effects of the presence of a nonlifting model are presented showing the resultant 
velocity component at different Mach numbers in Figure 5.6. The effect of the 
model is negligible although slight differences axe observed for suction. Piecewise 
linear approximation for suction and blowing agrees with the data. Since all the 
experiments were conducted at zero angle of attack, the effect of a lifting model 
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remains unknown. But, it is expected that influence coefficient matrices obtained 
at zero angle of attack would give a fairly good estimate for the pressure corrections 
at small angles of incidence. 



Figure 5.6: Effect of a nonlifting model on influence coefficients evaluated from 
side-wall pressure distribution. 
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5.1.4 Superposition 

The velocity change produced at any control point due to the pressure changes in 
more than one plenum was assumed to be equivalent to the sum of the velocity 
changes produced at that point by the same pressure changes in each plenum, if 
applied separately. A check of this assumption of superposition is given in Fig- 
ure 5.7. In this case, blowing is applied to plenum 4 and 6, first separately and 
then together. The velocity distribution for the combined plena, which was calcu- 
lated from side-wall pressure measurements, is compared to the superposed result. 
There is reasonable agreement although there are differences greater than 1 m/s at 
some points. These differences are attributed to the slight differences in freestream 
conditions and interpolation errors. 



X/C 


Figure 5.7: Superposition check for two plenum compartments for resultant 
velocity. 
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In Figure 5.8a , another superposition check using resultant velocity is presented. 
In this case, top and bottom plenum 6 are treated separately. Suction and blow- 
ing were first applied into the top plenum compartment. The same procedure was 
repeated for the bottom plenum while the others were kept closed. Influence coeffi- 
cients at control point 12 were calculated for each case and superposed to compare 
with those obtained when top and bottom plenum 6 were controlled together. 

Additional tests were performed to investigate the superposition in the presence 
of the model and at a higher Mach number. Influence coefficients at two control 
points upstream and downstream of plenum 6 were measured, and superposed using 
the same procedure outlined in the previous paragraph. Figures 5.8a and 5.8b 
illustrate that superposition for resultant velocity at a downstream control point is 
within the limits of the experimental error. Similar results were obtained at the 
upstream points as shown in Figure 5.8c. 



Figure 5.8a: Superposition check for one plenum compartment tor resultant 
velocity at M=0.5. 
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PLENUM NO. 6 M = 0.7 

o TP + BP 6 EXPERIMENT 
* TP 6 EXPERIMENT 
^ BP 6 EXPERIMENT 
o TP + BP 6 SUPERPOSED 
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Figure 5.8b: Superposition check for one plenum compartment for resultant 
velocity at M=0.7. 

CONTROL POINT NO. 6 



-.04 -.02 0 .02 .04 .06 .08 
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Figure 5.8c: Superposition check at an upstream point for one plenum compart- 
ment at M=0.7. 
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5.1.5 Linearity 

Over the Mach number range of interest, the influence coefficients are well approx- 
imated as linear for blowing at high Mach numbers. For suction, nonlinearity is 
clearly observed at low Mach numbers (Figures 5.5b, 5.6 and 5.8a ). As the Mach 
number increases, the linear approximation gets better for both suction and blowing 
(Figures 5.5a, 5.5b, and 5.6). In the present study, elements of the influence coeffi- 
cient matrices were approximated as piecewise linear for both suction and blowing. 
Since suction and blowing created different effects, two separate sets of influence 
coefficients were used in the construction of the influence coefficient matrices. In 
the present study, the effects of the nonlinear coefficients were not investigated the- 
oretically or experimentally since the linear approximation was found satisfactory 
for practical applications (See Section 5.2). 

5.1.6 Summary 

In the application of the convergence schemes, influence coefficient matrices were 
formed using the equations given in Chapter 2. Influence coefficients were approx- 
imated as piecewise linear in the present Mach number range. Complete sets of 
influence coefficient matrices were formed without the model and at M = 0.5 from 
the side-wall pressure measurements. In the case of LDV measurements, each set 
of matrices was constructed from influence coefficients obtained in the presence of 
the nonlifting model and at M = 0.5. The best linear fit was found for both suction 
and blowing at each control point. 

The influence coefficients obtained from side-wall pressure distribution are com- 
parable with those from streamwise velocity measurements using LDV. This sug- 
gests the possibility of using the same strategy for the unconfined flow calculations. 
Influence coefficient matrices determined for each flow variable are given in Ap- 
pendix C. 
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5.2 Application of the Convergence Schemes 

Adaptive-wall experiments were performed at Mach numbers ranging from 0.5 to 
0.75 and at angles of attack of 0°, 2°, and 4°. A brief summary of the test conditions 
is given in Table 2. 

The tunnel running time available for iterations was limited by the ejector which 
could operate continuously for only about 40 minutes at the maximum setting. At 
the end of each cycle, the wind tunnel was shut down so that the storage tank which 
drives the ejector, could be pressurized. The existing auxiliary system provided 
insufficient suction and blowing of air for complete control of the flow at high Mach 
numbers. 

Experiments with the LD V system were conducted at freestream Mach numbers 
up to M = 0.7; data quality at higher speeds was affected by condensation effects. 

In adaptive-wall experiments with the ventilated-wall tunnel, the test section 
Mach number varies when suction and blowing are applied in the plenum com- 
partments. During the experimental iterative process, after each iteration, the test 
section Mach number was adjusted back approximately to its initial value. In the 
one-step method, a first-order correction, deduced from calibration tests, was ap- 
plied to simulate this change in the test section Mach number. 

In this section, results obtained using the one- and two-level compatibility as- 
sessment methods with the iterative and the one-step convergence schemes will be 
presented in three parts. In the first part, typical results are given to establish 
the feasibility and the validity of the concept lasing the resultant velocity distribu- 
tion as a flow variable, as well as to show the applicability of the one-step method 
with the approximate influence coefficients in a supercritical lifting case. Secondly, 
the results obtained in a subcritical nonlifting case are presented to compare the 
one-level and two-level compatibility assessments along with the comparison of the 
iterative and the one-step methods using the streamwise and the resultant velocity 
distributions. In the last section, typical results are given to show the applicabil- 
ity and the effectiveness of the influence coefficients in subcritical and supercritical 
cases at various Mach numbers and model angles of attack. 
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As mentioned in Chapter 1, the finite number of plenum compartments is one of 
the major factors influencing the accuracy of adaptive-wall control. In the present 
study, a maxim um eight out of ten plenum compartments were used with the ex- 
ceptions being the first and the last pairs. Pressure measured at a location close 
to the middle section in each plenum compartment was assumed to be the plenum 
pressure in that compartment. In some of the examples, pressure distributions in 
the top and the bottom plenum compartments are presented for compaxison before 
and after the corrections are applied. In all the applications, there was initially no 
suction and blowing applied in the plena. After the final iteration, the status of 
the pressure control is illustrated for each compartment to indicate whether suction 
or blowing was applied. The size of the plenum compartments is drawn to scale, 
so that the influence regions can be easily seen. Since any pressure change in one 
plenum compartment induces changes in the others, pressure changes in plenum 
compartments 1 and 10 are due to the corrections applied in the rest of them. 


Applications with side-wall pressure measurements 

M 

0.5 

0.55 

0.6 

0.65 

a deg 

0, 2, 4 

0 

0, 4 

0, 2 

Scheme 

ite.; one-step 

iterative 


ite.; one-step 


M 

0.7 

0.72 

0.74 

0.75 

a deg 

0, 2, 4 

0 

0, 2 

0 

Scheme 

ite.; one-step 

ite.; one-step 

iterative 

iterative 


Applications with LDV 

M 

0.5 

0.65 

0.70 

a deg 

0, 2 

2 

0, 2 

Scheme 

ite.; one-step 

ite.; one-step 

ite.; one-step 

Flow 

Variable 

normal, and , 
streamwise vel. 

streamwise 

velocity 

normal, and , 
streamwise vel. 


Table 2: Test conditions in the experiments for the application of the conver- 
gence methods. 














































CHAPTER 5. EXPERIMENTS 


63 


5.2.1 Applications in the Supercritical Lifting Case 

Since the adaptive- wall concept has been developed primarily for transonic testing, 
the approach introduced in this research was first exploited in experiments where 
the model had a supercritical flow region on its surface. In some cases, nearsonic 
or transonic flow extended to the control levels before the wall boundary conditions 
were adjusted using the convergence schemes. 

5. 2. 1.1 Two-Level Iterative Scheme with Streamwise Velocity Distribu- 
tion 

The objective of this particular case was to use the streamwise velocity component 
as the flow variable in the two-level compatibility assessment and to establish a base 
for schemes with resultant velocity distribution. 

Since normal and streamwise velocity components were measured together, re- 
sultant velocity distributions at the control levels could be calculated and exper- 
imental steps could be simulated numerically by imposing the resultant “distur- 
bance” velocity distribution, u ' r , as the flow variable. The supercritical lifting case 
was chosen to find out the error introduced by this assumption, especially in the 
presence of high velocity gradients, in the flow field about the model. 

Initially, while there was no suction and blowing induced by the auxiliary system, 
velocity components were measured at M = 0.703 and a = 2.°. The two-level 
iterative method was applied with the streamwise velocity component only. In 
Figures 5.9a and 5.10a, experimental and calculated streamwise disturbance velocity 
distributions are shown at the upper and the lower levels respectively. Evidently, a 
shock wave penetrates the control levels in the upper half. However, in the lower 
half, the flow is subcritical and the effect of wall interference is much smaller. After 
the third iteration, wall interference was reduced substantially in the upper half, 
achieving a 60 % reduction in the rms error as shown in Figure 5.9b. Below the 
model, an almost perfect matching of the experimental and unconfined velocity 
distributions at level 2 was obtained as illustrated in Figure 5.10b. 
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UPPER LEVELS 
a=2. DEG. 



f ITERATION'S, M=0.698, RMS.=1.665 


□ MEASURED AT LEVEL 1 



b. 

| I I I 

- 2.0 - 1.0 0.0 1.0 2.0 

x/c 

Figure 5.9: Streamwise velocity distribution at the upper levels in the supercrit- 
ical lifting case, two- level iterative method. 
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ITERATIONS, M=0.703, RMS.=1.74 



□ MEASURED AT LEVEL 1 
o MEASURED AT LEVEL 2 
A UNCONFINED SOLUTION AT LEVEL 2 




ITERAT10N=3, M=0.698, RMS.=1.227 
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Figure 5.10: Streamwise velocity distribution at the lower levels in the super- 
critical lifting case, two-level iterative method. 
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The plenum pressure distributions were substantially different after three iter- 
ations compared to the initial case in which no suction and blowing was applied. 
Figure 5.11 shows the initial and the final pressure distributions of the each plenum. 
As expected, large pressure changes occurred in the plenum compartments in the 
vicinity of the model where the matching of the velocity distributions shows a large 
discrepancy, initially (Figure 5.9). 

M-0.700 a=2. DEG. 

o ITERATION: 0 

oo 



o 



x/c 


Figure 5.11: Plenum pressure distributions before and after the corrections ap- 
plied. 
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The initial and final measured model pressure distributions were compared with 
the Calspan data^ 45 ^ in Figure 5.12a. The effect of the wind tunnel wall interfer- 
ence is clearly observed on the model pressure distribution as the higher pressure 
coefficient and movement of the shock location between the initial and the final 
adapted cases. The Calspan data was taken in the Calspan eight-foot tunnel at a 
slightly higher Mach number compared to the present case, but it still provides a 
good comparison for our purpose. A transition strip used at the ten-percent chord 
location is believed to have been the cause for the suppression of the Calspan model 
pressure distribution. The Calspan data shows a local separation at ten-percent 
of the chord from the leading edge, while the present data also indicates a possi- 
ble local separation at x/C ~ 0.20. The present experiments were conducted at a 
nominal chord Reynolds number of 10 6 like the Calspan data. 


in 



Figure 5.12a: Comparison of the initial and final model pressure distributions 
with the Calspan data^ 45 b 
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Comparison of the initial and final model pressure distributions with the numer- 
ical solution from TAIR code* 47 * shows that the shock position is in good agreement 
after the corrections (Figure 12b). However, the numerical solution predicts a higher 
pressure coefficient at the leading edge compared to the experimental measurements. 
This is possibly due to the inviscid formulation used in the numerical code. 
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Figure 5.12b: Comparison of the initial and final model pressure distributions 
with the TAIR solution* 47 *. 
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Since the normal velocity component was measured along with the streamwise 
velocity, the changes in the normal disturbance velocity distributions could be com- 
pared. Figure 5.13 shows the initial and the final velocity distributions at the upper 
levels. It is observed that there is a considerable reduction in the rms error. 


UPPER LEVELS 

a=2. DEG. 




x/c 

Figure 5.13: Variation in the normal velocity distributions in the supercritical 
lifting case, two-level iterative method. 
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The experimental iterations outlined above were reconstructed numerically with 
the resultant velocity distributions which were obtained combining measured stream- 
wise and normal velocity components. The maximum error was found to be about 
1 % when the resultant “disturbance” velocity was assumed to be the streamwise 
velocity at the control points. Results show that rms errors are slightly lower, about 
4 %, compared to the scheme with the streamwise velocity distribution only. Oth- 
erwise, no significant differences were observed. The resultant velocity distributions 
are illustrated in Figures 5.14 and 5.15. By this simulation, it then becomes possible 
to establish the criteria with which the schemes can utilize the side-wall pressure 
distributions. 


UPPER LEVELS 
a= 2. DEG. 

ITERATIONS, M=0.703, RMS.=4.12 



x/c 


Figure 5.14a: Initial resultant velocity distributions computed from streamwise 
and normal velocity components at the upper levels. 
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Figure 5.14b: Resultant velocity distributions at the upper levels after the cor- 
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Figure 5.15: Initial and final resultant velocity distributions at the lower levels. 
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5. 2. 1.2 Two Level Iterative Scheme with Side-Wall Pressure Distribution 

Velocity components were measured with the LDV system along the centerline of the 
test section away from the side-wall boundary layers. However, side-wall pressure 
measurements were subject to the viscous effects caused by the presence of boundary 
layers. In the previous example, it was demonstrated that the convergence of the 
adaptive-wall concept to the interference-free flow could be achieved if the iterations 
were carried out with the resultant velocity distribution. The remaining question 
is whether resultant velocity distribution may be estimated from side- wall pressure 
measurements. In the subcritical case, comparison of resultant velocity distributions 
from side-wall pressure measurements with those from LDV measurements, showed 
that the boundary layer assumption holds fairly well. At supercritical speeds or in 
a lifting case on the model, a complex flow field may exist on the side- walls caused 
by the presence of a shock wave. In that case, the side-wall pressure measurements 
may not give a good estimate for the velocity distribution. To explore the effects 
of supercritical flow and to validate the feasibility of the velocity simplification, the 
iterative method was applied with the side-wall pressure distribution. 

Figures 5.16a and 5.17a illustrate the initial velocity distributions at the upper 
and the lower levels, respectively. Measurements were taken at M = 0.702 and at 
a = 2.° Initially, a weak shock wave was present at the first level in the upper half. 
The comparison of this case to the one with the streamwise velocity distributions 
at the same Mach number and angle of attack (Figures 5.9a and 5.10a), shows 
that side-wall pressure measurements underestimate the velocity distributions at 
the control levels. Using the estimated resultant velocity distributions, the iterative 
scheme was applied and wall interference was almost eliminated at the upper levels 
after five iterations as shown in Figure 5.16b. About 70 % reduction in the rms 
error was achieved in the velocity distributions in the upper half. At the lower 
levels, some improvement is obtained in matching the unconfined flow solution to 
the measured velocity distribution over the model (Figure 5.17b). 




Figure 5.16: Resultant velocity distributions from side- wall pressure measure- 
ments at the upper levels in the supercritical lifting case, two-level iterative method. 
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a=2. DEG. LOWER LEVELS 




X/C 

Figure 5.17: Resultant velocity distributions from side- wall pressure measure- 
ments at the lower levels in the supercritical lifting case, two-level iterative method. 


Top and bottom plenum pressure distributions are given in Figure 5.18. Com- 
parison with the previous case (Figure 5.11) shows similar pressure adjustments in 
the plenum compartments directly above and below the model while the last ones 
differ slightly. As pointed out previously, the small changes in the downstream 
compartments do not affect the upstream disturbance velocity distribution. 
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M=0.700 a=2. DEG. 

ITERATION: 0 


o 



o 



x/c 

Figure 5.18: Plenum pressure distributions before and after the corrections. 

Model pressure distributions measured after the initial and the final iterations, 
are compared with the Calspan data^ 45 ) and the numerical solution from the TAIR 
code^ 47 ) in Figure 5.19. The initial position of the shock wave is at the 50 %-chord 
of the model. In the final iteration, it was moved to a position in good agreement 
with the other results. 


CHAPTER 5. EXPERIMENTS 


76 


a— 2. DEG. 


ITERATIONS, M=0.702 : 




Figure 5.19: Comparison of the initial and final model pressure distributions 
with the Calspan data^ 45 ' and TAIR solution^ 47 ). 
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5. 2. 1.3 Two-Level One-Step Scheme with Side-Wall Pressure Distribu- 
tion 

After having shown the validity of using the side-wall pressure distribution in the 
adaptive-wall concept, the next step was to establish the applicability of the one- 
step scheme using the linearized influence coefficients. 

Figures 5.20a and 5.21a show the initial resultant velocity distributions which 
were obtained from side-wall pressure measurements at M = 0.699 and a = 2° 
at the upper and the lower levels respectively. Cumulative pressure corrections 
were determined by the one-step convergence scheme and they were applied in the 
plenum compartments as the final experimental step. 


UPPER LEVELS 
a =2. DEG. 



x/c 


Figure 5.20a: Initial resultant velocity distributions from side- wall pressure mea- 
surements at the upper levels in the supercritical case, two-level one-step method. 
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Figure 5.20b shows the velocity distributions at the upper levels after one itera- 
tion. Although the rms error decreased considerably (about 58 %), a small amount 
of interference remained in the velocity distribution at a portion of the control 
level over the quarter-chord point of the model. However, a small variation is also 
observed at the corresponding lower level points (Figure 5.21b). 


UPPER LEVELS 
a=2. DEG. 


ITERATIONS, M=0.70, RMS.=1.49 
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Figure 5.20b: Resultant velocity distributions from side-wall pressure measure- 
ments at the upper levels in the supercritical case after the corrections, two-level 
one-step method. 
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LOWER LEVELS 
a=2. DEG. 
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Figure 5.21: Resultant velocity distributions from side-wall pressure measure- 
ments at the lower levels in the supercritical case before and after the corrections 
using the one-step convergence scheme. 

Comparison of the top and bottom plenum pressure distributions after the first 
iteration (Figure 5.22) with those from the application of the iterative scheme (Fig- 
ure 5.18) reveals that plenum pressures were altered in the right direction. However, 
it also indicates that further adjustments are required in some plenum compart- 
ments. 
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Figure 5.22: Plenum pressure distributions before and after the corrections using 
the one-step convergence scheme. 

The discrepancy observed in the comparison of the measured model pressure 
distribution with the experimental 45 ) and the computational data (47 \ (Figure 5.23) 
is due to the residual interference that remained in the velocity distribution because 
of incomplete pressure adjustments. The position of the shock is slightly off toward 
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the trailing edge and the pressure coefficient is higher at x/C ~ 0.10 compared to 


the well-converged case in Figure 5.19. 



m 



Figure 5.23: Comparison of the initial and final model pressure distributions 
with the Calspan data* 45 * and TAIR solution* 47 *. 
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5. 2. 1.4 Discussion 

These results suggest that the estimated resultant velocity distribution from side- 
wall pressure measurements and the approximate calculation of unconfined flow in 
the schemes with side- wall pressure distribution may be acceptable for adaptive- wall 
experiments. 

The comparison of the iterative scheme (Figures 5.16 and 5.17) and the one- 
step scheme (Figures 5.20 and 5.21) reveals that the number of iterations could be 
reduced at least by three in this particular case with the application of the one- 
step convergence scheme. In this supercritical lifting case, it was also shown that 
convergence can be achieved with the use of linearized influence coefficients. 

As shown in Figure 5.24a, in the presence of a shock wave that has penetrated the 
control level, the unconfined solution using the transonic small disturbance equation 
predicts a different velocity distribution at level 2 compared to the solution obtained 
from the linear equation. 



x/c 


Figure 5.24a: Comparison of the linear and nonlinear unconfined flow solutions 
at level 2 with the measured data before the corrections in the supercritical case. 
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After the final iteration, linear and nonlinear solutions are in agreement since 
there is no shock wave present (Figure 5.24b). In spite of the initial presence of a 
shock wave at the control levels, the linearized equations were used to calculate the 
unconfined flow field and iterations were carried out successfully. 



x/c 


Figure 5.24b: Comparison of the linear and nonlinear unconfined flow solutions 
after the corrections. 

Approximate calculations of the unconfined velocity distribution in the nonlin- 
ear flow field and plenum pressure corrections with the use of linearized influence 
coefficients are probably the major reasons for not obtaining a complete conver- 
gence in one-step. It is inferred that one more iteration would have achieved a 
better convergence to the interference-free flow as obtained in the iterative case. 

5.2.2 Applications in the Subcritical Nonlifting Case 

The purpose of the applications in this section is to demonstrate that the two-level 
iterative and one-step convergence schemes can easily be implemented with the 
simplified flow variable. With the advantage of having fast and reduced number 
of iterations, similar results can be obtained as compared to those of the one-level 
two-component iterative scheme. In order to compare the results properly, all the 
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experiments were carried out at M ~ 0.7 and a = 0.° 

5. 2. 2.1 One-Level Iterative Scheme with Disturbance Velocity Compo- 
nents 

As stated previously the one-level two-component compatibility assessment was in- 
troduced by Sears* 15 ), and verified experimentally in several studies* 8,18 ). In this 
example, velocity measurements were taken at level 1. Figure 5.25a shows the ini- 
tial measured and calculated velocity distributions. Experiments were stopped at 
the fifth iteration because of time limitations and since little gain was apparent 
after the fourth iteration. Some discrepancies are observed between the unconfined 
solution and the measured streamwise velocity distribution after the final iteration 
as illustrated in Figure 5.25b. Some of the error is probably due to the slight varia- 
tion in the freestream velocity and the sensitivity of the integration to the changes 
in the normal velocity component. The measured model pressure distribution is il- 
lustrated in Figure 5.26. After the final iteration, it is generally in good agreement 



Figure 5.25a: Initial streamwise and normal disturbance velocity distributions 
in the subcritical nonlifting case, one-level iterative method. 
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Figure 5.25b: Streamwise and normal disturbance velocity distributions in the 
subcritical nonlifting case after the corrections, one-level iterative method. 



Figure 5.26: Comparison of the initial and final model pressure distributions 
with an experimental data base^ 46 *. 
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5. 2. 2. 2 Two-Level One-Step and Iterative Schemes with Streamwise Ve- 
locity Distribution 

In the application of the iterative scheme, initial velocity distributions were taken at 
M = 0.702 and zero angle of attack (Figure 5.27a). A fairly good convergence to the 
interference-free flow was achieved after 2 iterations as illustrated in Figure 5.27b. In 
the second iteration, 10 control points were used for the wall interference assessment 
at level 2 while the velocity measurements were made at 15 points at level 1 to 
eliminate the tail-effects as much as possible for the calculation of the unconfined 
flow. After the second iteration, the rms error decreased about 58 % and the 
disturbance velocities at level 1 and level 2 were reduced. Two more iterations were 
carried out to improve the rms error, but the gain was minimal and the velocity 
distributions and the model pressure coefficient remained approximately the same. 
This is probably due to inexact application of the small pressure corrections. 

UPPER LEVELS 

a= 0. DEG. 
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Figure 5.27a: Initial streamwise velocity distributions in the subcritical nonlift- 
ing case, two-level iterative method. 
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UPPER LEVELS 
a= 0. DEG. 
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Figure 5.27b: Streamwise velocity distributions in the subcritical nonlifting case 
after the corrections, two-level iterative method. 


In Figure 5.28, model surface pressure distributions after the zeroth and final 
iterations are compared with the experimental data taken at ONERA^ 46 ^. The effect 
of the wall interference is clearly observed between the initial and the final iterations. 
The converged case after the second iteration agrees well with the data base since 
the interference was eliminated rather well compared to the previous case. 



Figure 5.28: Comparison of the initial and final model pressure distributions 
with the ONERA data^ 46 K 

The one-step convergence scheme was applied at the same flow conditions, that 
is, at about the same Mach number and the angle of attack. The cumulative pressure 
corrections were determined using the one-step algorithm. After the first iteration, 
the test section Mach number was adjusted to 0.702 from 0.744 since it varied when 
suction and blowing were applied. Figure 5.29 illustrates the velocity distributions 
before and after the pressure corrections. 


UPPER LEVELS 
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The measured model pressure distribution is compared with the ONERA data 
in Figure 5.30 and a good agreement is observed although a perfect matching of 
the velocity distributions at level 2 could not be obtained in one iteration. In these 
last two cases, the discrepancies between the unconfined solution and the measured 
velocities at level 2 are observed above the quarter-chord point of the model which is 
located at X/c = 0 (Figures 5.27 and 5.29). One reason is probably the incomplete 
control caused by the finite number of plenum compartments over this region with 
the high velocity gradient. 



Figure 5.30: Comparison of the initial and final model pressure distributions 
with the ONERA data^ 46 *. 
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5. 2. 2. 3 Two-Level One-Step and Iterative Schemes with Side- Wall Pres- 
sure Distribution 

The one-step and iterative convergence schemes were applied with the resultant 
velocity distributions determined from side-wall pressure measurements. 

Figure 5.31a illustrates the initial velocity distributions at the control levels at 
M = 0.706 and zero angle of attack for the case of the one-step method. After 
the application of cumulative pressure corrections which were determined using 
the one-step algorithm, wall interference was almost eliminated as illustrated in 
Figure 5.31b. Due to the suction and blowing applied in the plenum compartments, 
the test section Mach number increased to 0.715 and for comparison purposes with 
the iterative case, the model surface pressure distribution was taken at that Mach 
number. 

UPPER LEVELS 
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Figure 5.31a: Initial resultant velocity distributions from side- wall pressure mea- 
surements in the subcritical nonlifting case, two-level one-step method. 
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Figure 5.31b: Resultant velocity distributions from side-wall pressure measure- 


ments after the one-step iteration in the subcritical nonlifting case. 



Figure 5.32: Comparison of the initial and final model pressure distributions 
with the computational data^ 48 ). 
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The iterative scheme was started at M = 0.726 as shown in Figure 5.33a and the 
flow initially was nearsonic at the airfoil surface. As presented in Figure 5.33 a, b, c 
and d, convergence to the interference-free flow was achieved in three iterations. 
After two iterations, wall interference was eliminated significantly. The last iteration 
further reduced the rms error for a better convergence. Initial and final model 
pressure distributions were compared with the computational data compiled by 
Lock* 48 *. Figures 5.32 and 5.34 show that similar results were obtained after the 
final iterations for both methods. 


UPPER LEVELS 
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Figure 5.33a: Initial resultant velocity distribution from side-wall pressure mea- 
surements in the subcritical nonlifting case, two-level iterative method. 
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UPPER LEVELS 

r ITERATIONS, M=0.727, RMS.=1.204 a=0. DEG. 



ITERATION=2, M=0.725, RMS.=0.697 

o MEASURED AT LEVEL 1 



- 2.0 - 1.0 0.0 1.0 2.0 
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Figure 5.33b,c,d: Convergence sequence in the application of the iterative con- 
vergence scheme in the subcritical nonlifting case. 
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o a—0. DEG. 



Figure 5.34: Comparison of the initial and final model pressure distributions 
with the computational data* 48 *. 


5.2.3 Additional Applications with the Side- Wall Pressure 
Distribution 

5. 2. 3.1 Subcritical Nonlifting Case 

Because of the blockage effect, the maximum Mach number attained in the test 
section was approximately limited to M = 0.75 with a nonlifting model if there 
was no wall control. At these flow conditions, the iterative scheme was applied 
with the side- wall pressure distribution. Initially, the flow was supercritical at some 
points of both control levels. Although the model was not at incidence, wall effects 
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were very severe on the velocity distributions and the model pressure coefficients 
at this Mach number (Figure 5.35a and 5.36). Convergence to the interference free 
flow was achieved in three iterations (Figure 5.35b). The rms error in the velocity 
distributions dropped from 4.799 to a final value of 0.806. 




x/c 

Figure 5.35: Resultant velocity distributions from side-wall pressure measure- 


ments in the subcritical nonlifting case before and after the corrections, two-level 
iterative method. 
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The comparison of the measured model pressure coefficients with the experimen- 
tal data base* 46 ) before and after the convergence, shows the drastic improvement 
in Figure 5.36. 


a= 0. DEG. ITERATIONS, M=0.751 : 



Figure 5.36: Comparison of the initial and final model pressure distributions 
with an experimental data base* 46 ). 

5. 2. 3. 2 Subcritical Lifting Case 

As noted before, although the influence coefficients were evaluated at one flow con- 
dition, they were used in determining the pressure corrections at different Mach 
numbers and model incidences. At a = 2.° and M = 0.654, the effect of wall in- 
terference initially was less severe on the model pressure coefficient (Figure 5.38) 
compared to the lifting case at M = 0.7 (Figure 5.12). In this case, excellent results 
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were achieved in five iterations which were carried out with side-wall pressure mea- 
surements (Figures 5.37a and 5.37b). After the final iteration, the measured model 
pressure distribution agrees well with the data base^ 48 * as illustrated in Figure 5.38. 



x/c 

Figure 5.37: Initial and final resultant velocity distributions from side- wall pres- 
sure measurements in the subcritical lifting case, two-level iterative method. 
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a= 2. DEG. ITERATION=0, M=0.654 : 



x/C 


Figure 5.38: Comparison of the initial and final model pressure distributions 
with the computational data* 48 *. 

5. 2. 3. 3 Subcritical Lifting Case 

In one specific case, model incidence was set to the maximum a = 4.° and it 
was considered sufficiently high for the purpose of these experiments. It took six 
iterations to arrive at the unconfined flow conditions at M = 0.6 while the rms error 
in the velocity distribution was reduced by about 70 % (Figures 5.39a and 5.39b). 
The flow was supercritical on the airfoil surface. Figure 5.40 shows the initial and 
final measured model pressure distributions and the experimental data base* 46 * for 
comparison. In general, the agreement is good and discrepancies are attributed to 
the differences in the angle of attack. 



Figure 5.39: Initial and final resultant velocity distributions from side-wall pres- 
sure measurements in a subcritical lifting case, two-level iterative method. 
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o a=4. DEG. 



x/C 


Figure 5.40: Comparison of the initial and final model pressure distributions 
with the ONERA data* 46 ). 

5. 2. 3. 4 Supercritical Lifting Case 

In extreme cases, such as M = 0.7 and a = 4.°, wall effects were more pronounced 
and complete wall control could not be provided because of insufficient suction 
and/or blowing capability. However, to show the capabilities of the schemes with 
the side-wall pressure distribution in such a severe case, the iterative and the one- 
step schemes were implemented with the available suction and blowing. Initially, a 
great amount of wall interference was present as shown in Figure 5.41a. Although it 
was considerably reduced in five iterations, discrepancies in the velocity distribution 
over the quarter-chord of the model could not be eliminated because of the lack of 
the sufficient wall control, Figure 5.41b. There was little interference at the lower 
control levels and it was improved slightly as shown in Figure 5.42. 



Figure 5.41: Initial and final resultant velocity distributions from side-wall pres- 
sure measurements at the upper levels in a supercritical lifting case, two-level iter- 


ative method. 
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The same flow conditions were obtained with both schemes after the final it- 
eration as presented in Figures 5.43 and 5.44. The shock location on the model 
upper surface was moved considerably and the pressure distribution was improved 
on the lower surface in spite of the limited wall control and the use of the linear 
flow analysis. 



x/c 


Figure 5.42: Initial and final resultant velocity distributions from side-wall pres- 
sure measurements at the lower levels in a supercritical lifting case, two-level iter- 
ative method. 
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Figure 5.43: Comparison of the initial and final model pressure distributions 


with the ONERA data* 46 * _ iterative method. 



with the ONERA data* 46 ) _ one-step method. 



Chapter 6 

Comparison of the Convergence 
Schemes 


One of the primary objectives of this study is to compare the different convergence 
schemes applied with different flow parameters. Based on the results presented 
in the preceeding chapter and the available data from other sources, comparisons 
made in three different categories are given in the following sections. It should 
be noted that the present experiments are not designed to create a data base for 
the pressure distributions for a NACA 0012 airfoil section. But instead, model 
pressure distributions were measured and compared with the other experimental and 
numerical data to demonstrate the consistency of the results and the applicability 
of the methods presented in this study. Moreover, in spite of the abundance of 
the experimental data on NACA 0012 airfoil, discrepancies between different wind 
tunnel results are not uncommon and a careful analysis is required to ensure the 
accuracy and quality of the data^ 50 \ 


6.1 One-Level versus Two-Level Convergence Scheme 

First, results from the applications of the one-level method^ 15 ) and the two-level 
method^ 26 ) will be compared at similar flow conditions. Figure 6.1 illustrates the 
comparison of the measured model pressure distributions after the final corrections. 
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ONE-LEVEL METHOD : 



Figure 6.1: Comparison of the model pressure coefficients from the application 
of the one-level and the two-level iterative methods and with the ONERA data* 46 ). 


Both methods were applied iteratively in this nonlifting case, and similar pres- 
sure coefficients were obtained except at the regions between 10 % and 30 % of the 
chord. This discrepancy is partly due to the imperfect matching of the velocity dis- 
tributions at the control level (Figure 5.25b). The maximum difference between the 
two pressure coefficients is about 0.06 (relative maximum variation=12 %, average 
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distribution while both velocity components were required in the one-level method. 
In the application of the one-level method, the number of iterations is considerably 
higher compared to the other method. Both results were also compared with ex- 
perimental data^ 46 ) in which height/chord ratio is 3.71 and Re = 3.79 x 10 6 . In the 
present experiments, height /chord= 1.44 and Re ~ 10 6 . The pressure distribution 
from the two-level method is found to be in good agreement with the data base. 


6.2 Convergence Schemes with Streamwise Ve- 
locity versus those with Side- Wall Pressure 
Distribution 

Side-wall pressure distribution was used as the flow variable in the application of 
the adaptive- wall concept to simplify the instrumentation and to reduce the testing 
time. To demonstrate the feasibility of this approach, the iterative and the one- 
step convergence schemes were applied using both streamwise velocity and side- wall 
pressure distributions. The results presented in Chapter 5, have shown that the two- 
level compatibility assessment method can be successfully applied with the side-wall 
pressure distribution and converges to unconfined flow conditions. 

Supercritical Lifting Case, (M ~ 0.7, a = 2.°) 

In Figure 6.2, the comparison of the measured model pressure distributions shows 
the consistency of the results from both applications in the supercritical lifting case. 
Both experiments were conducted at about the same Mach number and angle of 
attack. Although the number of iterations is higher in the experiments with the 
side-wall pressure measurements, the testing time for each iteration was about 5 
and 10 times shorter, in nonlifting and lifting cases respectively, compared to the 
experiments with the LDV measurements (based on same number of control points 
including on-line analysis and adjustment of the plenum pressures). 
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Both pressure measurements were compared with the data from the Calspan 
8-foot tunnel (Re = 10 6 and transition strip at x/c=0.10 at the leading edge). 
Although it was taken at a slightly higher Mach number, it provides a good com- 
parison for our purpose. The shock location seems to be in good agreement and 
boundary layer trip is probably the cause of the local separation and the lower 
pressure coefficient at the leading edge. 


ITERATIVE METHOD 

a=2. DEG. USING STREAMWISE VELOCITY : 



Figure 6.2: Comparison of the model pressure coefficients obtained from the 
application of the iterative methods using the streamwise velocity and the side-wall 
pressure. 
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Another comparison is given in Figure 6.3 for the same results with other data 
taken in the NASA Langley 0.3 meter cryogenic adaptive wind tunnel^ 51 ). This 
tunnel has a flexible-wall test section and the presented data was taken at M=0.70 
and Re=3.02 10 6 . Model angle of attack was given as 1.976° and it has a boundary 
layer trip at x/c=0.05 . In this case, the model incidence is slightly smaller and 
shock location is nearer the leading edge compared to our data. However, top and 
bottom pressure coefficients are found to be in good overall agreement. 


a =2. DEG. 


ITERATIVE METHOD 



Figure 6.3: Comparison of the model pressure coefficients obtained from the 
application of the iterative methods using the streamwise velocity and the side-wall 
pressure, and with the LANGLEY flexible-wall adaptive wall data^ 51 h 
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Subcritical Lifting Case, (M ~ 0.65, a = 2.°) 

Model pressure distributions taken at similar Mach numbers (M=0.649 for the ex- 
periments with streamwise velocity distribution and M=0.636 for the side- wall pres- 
sure measurements) are compared in Figure 6.4 . It can be seen that the pressure 
measurements for the lower surface are in good agreement, but the data on the 
upper surface differ, especially at the leading edge (the maximum difference=0.136, 
relative maximum variation~ 14 %, average variation on the upper surface^ 7 %). 



Figure 6.4: Comparison of the model pressure coefficients obtained from the 
application of the iterative methods using the streamwise velocity and the side-wall 
pressure and with the computational data^ 48 h 
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The experiments with the side-wall pressure distribution showed a perfect match- 
ing in the velocity distributions as shown in Figure 5.37 . A good convergence was 
also achieved in the application of the iterative scheme using the streamwise ve- 
locity. The discrepancies may be due to the slight differences in the setting of the 
model incidence. The present data are compared with the exact solution by Sells as 
reported by Locld 48 ) in Figure 6.4 and with the experimental data from the Calspan 
8-foot tunnel^ 45 ) in Figure 6.5 . The discrepancies are attributed to the differences 
in the Mach number, and the angle of attack. 


in 


a-Z. DEG. 


ITERATIVE METHOD 



x/C 


Figure 6.5: Comparison of the model pressure coefficients obtained from the 
application of the iterative methods using the streamwise velocity and the side-wall 
pressure and with the CALSPAN data( 45 b 
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6.3 Iterative Scheme versus One-Step Scheme 

The purpose of the one-step method is to obtain the unconfined flow conditions 
in a single iteration if possible, or otherwise to minimize the number of iterations. 
Experiments were conducted to demonstrate the relative merits of the iterative 
scheme and the one-step convergence scheme introduced in this study. 

Subcritical nonlifting case with streamwise velocity 

In Figure 6.6, the comparison of the measured model pressure distributions from the 
experimental iterative and one-step convergence schemes shows that similar results 
were achieved in one iteration. 



Figure 6.6: Comparison of the model pressure coefficients obtained from the 
application of the iterative and the one-step methods using the streamwise velocity. 
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Experiments were carried out at M ~ 0.7 and at zero angle of attack. Aside 
from the slight variations caused by the differences in the Mach numbers, the model 
pressure coefficients after the final iterations are found to be in good agreement. The 
present pressure data is also compared with the experimental data from ONERA^ 46 ) 
in Figure 6.6 and with the data taken in the NASA Langley 8-foot transonic pressure 
tunnel* 49 ^ in Figure 6.7. The ONERA data which was discussed previously is usually 
in good agreement with the present results in spite of the small discrepancies. 



Figure 6.7: Comparison of the model pressure coefficients obtained from the 
application of the iterative and the one-step methods using the streamwise velocity. 
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The variation between the upper and the lower surface pressure coefficients at 
0.10 < x/c < 0.30 in Harris’ data* 49 ) in Figure 6.7 would probably stem from the 
differences in the angle of attack (Re~ 3.0 x 10 6 in both data base). Otherwise, the 
present pressure data agrees reasonably well with Harris’ data over the rest of the 
chord length. 

Subcritical Nonlifting Case with Side-Wall Pressure Distribution 

In this case, the iterative and the one-step convergence schemes were applied using 
the side-wall pressure distribution at slightly higher Mach numbers compared to 
the previous example. Experiments with both schemes produced similar pressure 
coefficients as illustrated in Figure 6.8 . As seen from the comparison, convergence 
to the unconfined flow was achieved in one iteration using the one-step method 
while three iterative steps were required to arrive at the same conditions using the 
iterative convergence scheme. The small discrepancies could possibly stem from 
the differences in the Mach numbers. Experimental data was compared with the 
numerical solution reported by Lock* 48 ). The comparison is reasonably good over 
the entire chord length. Lock’s data in general is in good agreement with the 
Calspan data* 45 ) (M=0.725, a = 0.°) although it was not presented here. 

Supercritical Lifting Case with Side- Wall Pressure Distribution 

Experiments using the one-step convergence scheme did not produce the same re- 
sults in one iteration, however the effect of the wall interference on the model 
pressure distribution was greatly reduced as seen in Figure 6.9 . The pressure coef- 
ficient is slightly higher in the region 0.025 < x/c < 0.30 and the shock position is 
about 5 % further downstream compared to the iterative case. Imperfect matching 
in the velocity distributions (Figure 5.20b) in the same region is probably the main 
reason for the discrepancy. One more iteration would yield similar results to those 
obtained using the iterative scheme. Although unconfined flow conditions could 
not be achieved in one iteration, the number of iterations was reduced drastically 
compared to the iterative method. In Figure 6.9, the Calspan data is included for 
comparison with the measured model pressure distributions. 
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Figure 6.8: Comparison of the model pressure coefficients obtained from the 
application of the iterative and the one-step methods using the side- wall pressure. 

Supercritical Lifting Case with Side- Wall Pressure Distribution 

At M=0.7 and a = 4°, wall interference effects were severe and the available wall 
control was not sufficient to provide the necessary adjustments. However, the it- 
erative and the one-step convergence schemes were applied with the available wall 
control. In Figure 6.10, the comparison of the model pressure coefficients demon- 
strates the ability of the one-step method to produce similar corrections in fewer 
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iterations compared to the iterative scheme. Results are also compared with exper- 
imental data^ 4 " . Except for the shock locations, the present data is in agreement 
with the experimental data from ONERA. The discrepancy is obviously dependent 
on the incomplete wall corrections in the present experiments. 


a=2. DEG. ITERATIVE METHOD : 


in 



Figure 6.9: Comparison of the model pressure coefficients obtained from the 
application of the iterative and the one-step methods using the side- wall pressure. 
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(X-A. DEG. 


ITERATIVE METHOD : 


in 



Figure 6.10: Comparison of the model pressure coefficients obtained from the 
application of the iterative and the one-step methods using the side-wall pressure. 


6.4 Testing Time Comparisons 

In this section, a comparison of the typical testing times in the application of the 
convergence schemes is given to determine the relative time efficiency of the pre- 
ferred technique, the one-step convergence scheme using the side-wall pressure. 
Testing time for each iteration depended mostly on the measurement technique. 
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Velocity measurements using the LDV system usually took about 20 minutes for a 
nonlifting case (a total of 30 control points at two control levels) and about 40 min- 
utes for a lifting case. For a lifting case, velocity measurements had to be taken in 
the upper and lower control levels since the flow was not symmetrical with respect 
to the model. Collection of the LDV data consisted of data taking, processing and 
moving the LDV system along the control levels. Side- wall pressure measurements 
and data processing took about 1 minute regardless of model configuration. Ap- 
plication of the pressure corrections in the plenum compartments contributed an 
average of 5 minutes while the calculation of unconfined flow and pressure correc- 
tions were negligible, about 0.1 minute. 

In a typical application using the iterative method, the total testing time after 
three iterations was about 75 minutes using the LDV system. The same application 
using the side- wall pressure measurements took about 18 minutes, reducing the total 
testing time by a factor of approximately 5. Application of the one-step convergence 
scheme took approximately 8 minutes using the side- wall pressure measurements, 
which reduced the total testing time by a factor of approximately 10. 



Chapter 7 

Conclusions and 
Recommendations 


7.1 Conclusions 

A fast and reliable method for the elimination of the wind tunnel wall interference 
effects has been demonstrated. The main conclusions that have been arrived at in 
connection with the use of this method are as follows: 

1. The experiments demonstrated that it is feasible to use the side- wall pressure 
distribution as a flow variable in the adaptive-wall convergence schemes. The mea- 
surement technique requires simple instrumentation and reduces total testing time 
considerably. The convergence schemes applied using the side-wall pressure yielded 
accurate and consistent results. Measurements of the model pressure distributions 
agree well with the other experimental and computational data. 

2. Existence of linear influence coefficients independent of Mach number and 
presence of the model permits a fast correction scheme. Influence coefficient ma- 
trices, which were constructed using the superposition principle, could be used to 
determine the pressure corrections at different Mach numbers, with and without 
lift. 

3. The use of the less complicated side-wall pressure instrumentation instead 
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of the LDV measurement of velocity was found to reduce the testing time by a 
factor of approximately five. The one-step method of adapting the wall conditions 
to satisfy interference-free flow, using either the streamwise velocity or the side-wall 
pressure, further reduces the testing time by a further factor of two to three when 
compared to the iterative scheme. 

4. The adaptive-wall corrections determined by using the linearized flow equa- 
tions provided sufficient control to achieve convergence to interference-free flow for 
both the subcritical and mildly supercritical flows. 


7.2 Recommendations 

The recommendations for the future regarding adaptive-wall experiments are as 
follows: 

1. The use of influence coefficients should be further investigated to determine 
the limits of their applicability in the case of highly supercritical flows for which 
shock waves may extend to the tunnel walls. 

2. Automated control of suction and blowing to the plenum compartments 
should be incorporated for production testing in order to obtain faster, more precise 
wind tunnel wall adaptation. 

3. The number of plenum compartments should be increased especially over the 
region near the model to improve flow control. 



Appendix A 

Numerical Evaluation of the 
Cauchy Integral 


This section shows the removal of the singularity in the numerical evaluation of the 
Cauchy integral. On a singular strip: 

v . singularity V]+l 

! . Ae — ' . 




Xj+1 


The Cauchy integral equation for the evaluation of the streamwise disturbance 
velocity is given as: 


u(x,,) = -i^ 


+0 ° v(G,ri)dZ 

-oo £ — X 


Define lj = Xj+ j — Xj and £ = x — Xj. 

Assume that the normal disturbance velocity along the singular strip varies 
linearly as a function of £, that is; 


v((,v) = a£ + b 

The coefficients a and b can be evaluated using the boundary conditions given 
at the ends of the strip. They are given as follows; 
at £ = x-j v = vj and at £ = Xj+i v = Vj+j 
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Then; 


a = 


v j + 1 ~ Vj 
x j+i x j 


b = Vj — 


v j+1 - Vj 


x i 


x j + 1 - x j 

Hence, the functional relation for the normal velocity can be expressed as: 




Substituting this equation into the Cauchy integral equation, the following ex- 
pression is obtained: 




u,+i — tq . . 

v i + — 7 (* - x i) 

‘j 




At the singular point, a limiting value can be found about a circle assuming 
its radius, e, goes to zero while £ approaches x. Then, the equation above can be 
written as: 


u( *' n) = { Vi + V - Ji ir 1{z - *>>) (sstfjT rh + j'i rh>) 

As a result, the velocity about the singularity on the singular strip can be 
determined using the following nonsingular equation: 




”, + ^7 — — (i - x,)) ln(- — -) + (fj+i - ”,) 

V l j ) x 



Appendix B 


Transonic Small Disturbance 
Equation 


In this section, a finite difference formulation for the solution of the transonic small 
disturbance equation using Murmon-Cole method will be given. The transonic 
small disturbance equation can be written in the following form by applying the 
Mur man- Cole scheme: 


(! Vi ,) — 0 


where 






^ oo 4*+! ,j 4i-l,j 

?oo x t— i,j 


and 


Vi,j — 


0 if A it j > 0 

1 if A itj < 0 


Dxx ^ 2 f (f> j+l,j 4i,j _ \ 

Ax 2 (®t+l,j — — X i<j x i,j ~ X i~1,j / 

Dyy _ 2 f 4i,j + l ~ 4>i ,j _ ~ 4i,j-l \ 

Ay 2 9l ’ 3 ( Vi,j+1 - \ Vi,j + 1 - Vi,j Vi,j - Vi,j - 1 / 
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Ai-ij = 1 - Ml - (7 + 1 ) 


The Murman-Cole method becomes: 


Ml fc - 4>i-u 

9oo x i,j ~ x t-2,j 


a _ . , A , _ Ki-K-li 

— l,i) \ Xi+lj — x i,j x i,j — x i-\,j . 

2 (K? ~ fiSj 


+rf- hj Al hj 
4 


(Xj t j X{-2,j ) \ x i,j x t-lJ x i-l,j x i-2,j J 

2 ( imu-tn 1 0 

(y«ii+i — V y<.i+i — y«.i y«.i — y*.j- 1 / 

This equation can be arranged as follows: 

S.tfjii + Si^ 1 4- 54^ = -So<^ + W - - Sa^S,- 

where 


Si = (i - pWK 


+ 


' J (xi+ij \Xi+iJ X ij XiJ Xi—ijJ 

2 / 1 > 
(x, t y X t _2,j) \^t,j 1 »jV 

— 1 2 -^’ 

y«,j+i y«,i-i \y»,j+i y»,j y*,j y».j- 1, 


S 2 (1 ( (lj+w _ (xy - 

A »_ 2 f 1 | 1 

’ J 1,7 ( x iJ — X,_2,j) \£i,j — l,j X i-l,j x i—2,j 


S3 — MUj A h,s 


( X i,j x i-2,j) \ x i-l,j Xi-2,j, 
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S 5 = 


2 


Vi,j + 1 ViJ-1 



The numerical equations derived can easily be solved by using the Gauss-Seidel 
line iteration technique. The solutions of the transonic small disturbance theory 
are valid for small perturbations and for freestream Mach numbers only slightly 
different from sonic. In the present study, the computational domain does not 
contain the airfoil and the freestream Mach numbers do not exceed M = 0.75, so 
that the inviscid approximation and the assumption of small disturbances hold well. 
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rjwwob 


PL3 

PL4 

PL5 

PL6 

PL7 

PL8 

PL9 


U 

Suction 

20.4 

24.8 








Blowing 

2 

Suction 

22.3 

43. 

5.1 

12.5 







Blowing 

3 

Suction 

26.3 

56.7 

6. 

22.2 







Blowing 

D 

Suction 

27.6 

52.8 

11.6 

35.7 

3.1 

5. 






Blowing 

5 

Suction 

29.6 

48.9 

17.1 

49.1 

4.8 

10. 






Blowing 

D 

Suction 

30. 

42.4 

16.3 

47.5 

11.2 

20. 

2.7 

5.9 





Blowing 

n 

Suction 

27.1 

42. 

15. 

45. 

13. 

26.1 

3.5 

10.5 





Blowing 

8 

Suction 

28. 

45. 

13.3 

41.7 

14.8 

24.9 

5.7 

21.6 





Blowing 

9 

Suction 

25. 

48. 

13.6 

40.4 

16.6 

25.2 

7.9 

29.3 

2.4 

8.5 




Blowing 

D 

Suction 

26.5 

48.5 

13.7 

40.8 

17.6 

24.7 

11.3 

30. 

5.7 

16.6 




Blowing 

1 1 

Suction 

26.2 

44.2 

14.5 

42.7 

19.4 

23.6 

15.3 

28.4 

10.9 

25. 

12.4 

14. 



Blowing 

12 

Suction 

26.4 

40.4 

14.7 

39.9 

17.5 

24. 

13.5 

25.1 

10.7 

27.7 

23.4 

23.8 

6.8 

16.8 


Blowing 

13 

Suction 

27. 

39. 

14.9 

36. 

16.5 

24.2 

9.7 

20.7 

8.8 

19.6 

28.3 

35. 

11.8 

19.8 


Blowing 

1 4 

Suction 

24.6 

40.1 

14.3 

36. 

17.4 

24.4 

8.5 

20.8 

8.2 

17.5 

26. 

36.7 

19. 

24.3 

6.7 

9.6 

Blowing 

1 5 

Suction 

26.7 

39.5 

13.1 

38.9 

16.3 

25. 

8.7 

21.4 

8.2 

18.3 

25.5 

32. 

20. 

35.3 

12 . 

9.8 

Blowing 




Table C.l: Influence coefficient matrices for suction and blowing from streamwise 
velocity distribution at level 1 in the presence of a nonlifting model at M — 0.5, 
[AV]/[AP] = l(m/s)/psi}. 
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PL3 

PL4 

PL5 

PL6 

PL7 

PL8 

PL9 


Suction 20.1 
Blowing 23.3 


Suction | 23.1 5.1 


Blowing 38* 


Suction | 26. 


Blowing 52.8 18. 


Suction 28.2 11.3 

Blowing 49.3 39. 


Suction 27. 



Blowing 44. 


Suction | 26.9 13 

42.4 


\rmm\ 




Blowing 


Suction 


Blowing 


Blowing! 40. 


Suction | 26.6 15. 


irroni 


Suction 24.7 
Blowing 3 


Suction | 24.5 



11.3 

39.7 

LA 

U) 

LA 


15.6 

55. 

4.8 

9.3 


14.7 

11.3 

1 . 

52. 

24.7 

5. 

14.1 

13.6 

2. 

45.7 

28.2 

8. 

13. 

16.6 

4.5 

37.5 

26.2 

22.2 

15. 

17.8 

7.4 

33.9 

22.5 

28.4 

14.8 

19. 

10.5 

33.7 

22.8 

30. 

16.9 

16.7 

14.8 

36.5 

21.1 

26.4 

15. 

15.2 

13.5 

39.9 

21. 

22.2 

15.3 

15.7 

11.4 

34.6 

18.9 

18. 

15.2 

17.4 

9.8 

32. 

17.3 

20.5 

13.1 

19. 

9.5 

28. 

18.6 

21.8 



Table C.2: Influence coefficient matrices for suction and blowing from streamwise 
velocity distribution at level 2 in the presence of a nonlifting model at M = 0.5, 
[AV]/[AP) = { (m/s)/ps,}. 
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I 


he EESM 


Suction 


Blowing 


Suction tg.2 



Blowing! 40.7 11.1 


3 Suction 21.2 8. 



Blowing 44.4 24. 


Suction 22.2 11.6 

Blowing 42.2 40. 


Suction 22.3 13.5 8.9 

Blowing 41.1 47.5 21. 


Suction 22.4 14.1 13.6 

Blowing! 40.1 45. 39. 












Suction 

22.2 

14.1 

15. 

4.4 

1.8 

Blowing 

40. 

43. 

43.6 

13.6 

6. 



8 

Suction 

22.2 

14.1 

16.1 

5.6 

BUI 



Blowing 

40. 

41. 

45. 

18.6 

HkB 



g 

Suction 

22.2 

14. 

16.7 

6.8 

4.3 




Blowing 

40. 

40. 

43.7 

22. 

13.7 



D 

Suction 

22.2 

14. 

17.1 

7.4 

5.9 

4.7 


Blowing 

40. 

40. 

42. 

23. 

19.5 

6.8 


1 1 

Suction 

22.2 

14. 

16.9 

7.9 

7.7 

9.9 


Blowing 

40. 

40. 

39.4 

21.5 

23.7 

17.6 


1 2 

Suction 

22.2 

14. 

16.5 

7.9 

8. 

14.9 

5.2 

Blowing 

40. 

40. 

38.4 

20.2 

22.1 

33.2 

7.6 

1 3 

Suction 

21.4 

14. 

16.6 

7.7 

7.8 

17.7 

9.9 

Blowing 

39.5 

40. 

38.1 

20.2 

20.3 

39.2 

17.7 

1 4 

Suction 

21.4 

14. 

16.6 

7.7 

7.7 

18.5 

15.2 

Blowing 

39.5 

40. 

38. 

20.2 

19.1 

35.3 

32.1 

1 5 

Suction 

21.4 

14. 

16.4 

7.6 

7.7 

18.4 

18.4 


Blowing 

39.5 

40. 

38. 

20.2 

19.1 

32.7 

38.8 



Table C.3: Influence coefficient matrices for suction and blowing from side- 
wall pressure distribution at level 1 without a model at M = 0.5, [AV]/[AP] = 
[(m/ 5 )/psi]. 
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QHjnaaisj 

jTTiimiirtfnnl 



Suction 15.7 
Blowing 24.7 


2 

Suction 

20.1 


Blowing 

44.7 


Suction 

21.2 


Blowing 

46.8 


Suction 22.2 
Blowing 42.2 


PL3 PL4 PL5 PL6 PL7 PL8 PL9 



3.8 

9.5 




5 

Suction 

22.3 

14.2 

8.9 


Blowing 

40.6 

51. 

19.5 

£ 

Suction 

22.4 

14.5 

14.6 

D 

Blowing 

39.8 

45. 

44. 



m 

Suction 

22.2 

14.3 

16.0 

4.4 

1.3 

mi 

Blowing 

40. 

43. 

49.6 

13.5 

4.4 

8 

Suction 

22.2 

14.2 

16.9 

6.2 

2.5 

MLJSllin 

40. 

41. 

48. 

21.3 

8. 

Q 

Suction 

22.2 

14.1 

17.2 

7.3 

4.3 

7 

Blowing 

40. 

40. 

44.7 

25.4 

14.2 

D 

Suction 

22.2 

13.9 

17.1 

7.9 

6.3 

Blowing 

40. 

40. 

41.7 

24.6 

21.5 

1 1 

Suction 

22.2 

13.9 

16.9 

8.3 

8.4 

Blowing 

40. 

40. 

39.4 

21.8 

25.9 

12 

Suction 

22.2 

14.1 

16.5 

8. 

8.3 

Blowing 

40. 

40. 

38.4 

20.2 

22.1 

1 3 

Suction 

21.4 

14. 

16.6 

7.7 

7.8 

Blowing 

39.5 

40. 

38.1 

20.2 

20.3 

14 

Suction 

21.4 

14. 

16.6 

7.7 

7.7 

Blowing 

39.5 

40. 

38. 

20.2 

19.1 

1 5 

Suction 

21.4 

14. 

16.4 

7.6 

7.7 


Blowing 

39.5 

40. 

38. 

20.2 

19.1 



Table C.4: Influence coefficient matrices for suction and blowing from side- wall 
presure velocity at level 2 without a model at M = 0.5, [AV]/[AP] = [( m/s)/psi }. 
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Control Points 


r 


Plenum Com.' 


PL 2 


PL3 


PL4 


PL5 


PL6 


PL7 


PL8 


PL9 


1 


Blowing 


Suction 


- 5.57 

- 10.24 


Suction 


Blowing 


- 3.52 

- 5.45 


- 1.9 

- 7.31 


Suction 


Blowing 


- 2.79 

- 7.77 


Suction 


Blowing 


- 2.12 

- 6.3 


• 1.75 

- 5.22 


Suction 


Blowing 


- 2.81 

- 7.63 


Suction 


Blowing 


- 4.86 

- 5 . 


- 2.14 

- 4.36 


Suction 


Blowing 


- 3.73 

- 5.16 


-2, 


37 

- 4.56 


Suction 


Blowing 


- 4.36 

- 4.68 


- 3 . 


51 

5.95 


Suction 


Blowing 


- 3.12 

- 3.28 


- 4 . 


5 

- 7.04 


1 0 


Suction 


Blowing 


1.29 

- 6.35 


- 4.26 

- 2.93 


1 1 


Suction 


Blowing 


- 2 . 


22 

- 2.61 


- 7.47 

- 9.96 


1 2 


Suction 


Blowing 


- 6.65 

- 5.12 


- 3 . 


68 

- 3.29 


1 3 


Suction 


Blowing 


-7 


.98 

- 6.36 


14 


Suction 


Blowing 


- 4 . 


18 
- 5.36 


1 5 


Suction 

Blowing 


- 4.99 
- 5.26 I 


Table C.5: Influence coefficient matrices for suction and blowing from normal 
velocity distribution at level 1 in the presence of a nonlifting model at M = 0.5, 
[AF]/[AP] = [( m/s)/psi ]. 
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Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


|;JMHiliL.I 




Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowin 


Suction 


Blowing 


Suction 


Blowing 


Suction 


Blowing 



Table C.6: Influence coefficient matrices for suction and blowing from normal 
velocity distribution at level 2 in the presence of a nonlifting model at M = 0.5, 
[AV]/[AP] = [(m/s)/psi]. 
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